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Abstract: 

Extending the work by Bukhvostov, Frolov, Lipatov and Kuraev (BFLK) on the renor- 
malization of quasipartonic operators we derive a complete set of two-particle renormal- 
ization group kernels that enter QCD evolution equations to twist-four accuracy. It is 
shown that the 2^2 evolution kernels which involve "non-partonic" components of 
field operators, and, most remarkably, also 2 — > 3 kernels do not require independent 
calculation and can be restored from the known results for quasipartonic operators using 
conformal symmetry and Lorentz transformations. The kernels are presented for the 
renormalization of light-ray operators built of chiral fields in a particular basis such that 
the conformal symmetry is manifest. The results can easily be recast in momentum 
space, in the form of evolution equations for generalized parton distributions. 
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1 Introduction 



Higher-twist effects generically correspond to corrections to hadronic observables that 
are suppressed by a power of the (large) momentum transfer or the heavy quark mass. 
One application where such effects are phenomenologically relevant are high precision 
studies of the total cross section of deep- inelastic lepton-hadron scattering (DIS). In this 
case the leading-twist calculations within the standard DGLAP formalism are advanced 
to the next-to-next-to-leading order. Also lattice calculations of the moments of parton 
distributions with an accuracy at a percent level are becoming feasible. In this situation 
taking into account twist-four corrections that are suppressed by a power of the photon 
virtuality Q 2 proves to be increasingly important for the analysis of modern data, see 
e.g. [D [21 El III El E] . Twist-three effects are actively discussed in the context of exclusive 
and semi-inclusive reactions, e.g. deeply- virtual Compton scattering [TIE] and diffractive 
electroproduction of vector mesons [9] , single spin asymmetry in various reactions [TUl 
ITTl IT2"j [T3l [HI [151 [16], etc. One can expect that with the increasing accuracy of the 
experimental data the twist-four effects will start playing a role here as well. Another 
large field of applications are the studies of higher-twist hadron distribution amplitudes 
(e.g. [ITJCTSldSlEn]) that provide one with the important input to the so-called light-cone 
sum rules [2T1I22]. 

The theoretical description of higher-twist corrections is based on the Wilson Op- 
erator Product Expansion (OPE) and involves contributions of a large number of local 
operators. The corresponding leading-order coefficient functions are usually easy to cal- 
culate (e.g. twist-four contributions to DIS are known since many years [221121]) but the 
operator renormalization for the operators of twist-four and higher has not been studied 
systematically. Up to now, twist-four anomalous dimensions are only known for a subset 
of four-quark operators [22] and for a few quark-gluon operators of lowest dimension 
(e.g. [26K2312HI 29]). In addition, the structure of the most singular parts of the mixing 
kernels for small values of the Bjorken variable that are relevant for the contribution of 
two-pomeron cuts in high-energy scattering processes was considered in [301 EI] • 

A general formalism was developed by Bukhvostov, Frolov, Lipatov and Kuraev 
(BFLK) [32] for the special class of so-called quasipartonic operators that are built 
of "plus" components of quark and gluon fields. For each twist, the set of quasipartonic 
operators is closed under renormalization and the renormalization group (RG) equation 
can be written in a Hamiltonian form that involves two-particle "interaction" kernels, 
cf. Fig. that can be expressed in terms of two-particle Casimir operators of the 
collinear subgroup SL(2, M) of the conformal group. In this formulation symmetries of 
the RG equations become explicit. Moreover, the corresponding three-particle quantum- 
mechanical problem turns out to be completely integrable for a few important cases, and 
in fact reduces to a Heisenberg spin chain ]33j. An almost complete understanding 
achieved at present of the renormalization of twist-three operators is due to all these 
formal developments, see [3H [35] for a review and further references. 

The goal of our study is to generalize the BFLK approach to the situation where not 
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Figure 1: Schematic structure of one-loop renormalization group kernels in QCD 



all contributing operators are quasipartonic, as it proves to be the case starting with 
twist four. On this way, there are two complications. 

First, the number of fields ("particles") is not conserved. To one-loop accuracy, the 
mixing matrix of operators with a given twist has a block-triangular structure as the 
operators with less fields can mix with ones containing more fields but not vice versa. 
Operators with the maximum possible number of fields are quasipartonic. 

Second, operators involving "minus" and "transverse" derivatives and/or field com- 
ponents must be included. The problem is that transverse derivatives generally do not 
have good transformation properties with respect to the SL(2, M) group. In concrete 
applications it may be possible to get rid of such operators using equations of motion 
(EOM) and exploiting specific structure of the matrix elements of interest, e.g. if there 
is no transverse momentum transfer between the initial and the final state. The main 
problem as far as the operator renormalization is concerned is that after this reduction 
the conformal symmetry becomes obscured. In the work [36J we have suggested a differ- 
ent, general approach based on the construction of a complete conformal operator basis 
for all twists. In this basis, the SL(2, M) symmetry of the RG equations is manifest. 

Another important observation was made in [37] where it was shown that the diagonal 
part of one-loop QCD RG equations (for arbitrary twist) can be written in a Hamiltonian 
form in terms of quadratic Casimir operators of the full conformal group SO(A, 2). This 
implies that the 2 — > 2 kernels of the type shown in Fig. [TJd can be obtained from the 
BFLK kernels in Fig. [T^l by the corresponding replacement. Our first goal will be to 
work out the corresponding relations using the conformal basis of Ref. for generic 
twist-four operators. 

The main new contribution of this work is the calculation of 16 independent 2^3 
kernels corresponding to operator mixing of one "partonic" and one "non-partonic" field 
in three-particle quasipartonic operators, shown schematically in Fig. [Ifc. To this end we 
suggest a new technique based on application of Lorentz transformations (translations 
and rotations in the transverse plane) to the "diagonal" 2^2 kernels, which bypasses 
calculation of Feynman diagrams. 

Combining the kernels in Fig. []Ji,b,c one obtains a complete set of building blocks 
for the renormalization of twist-4 operators that involve at most one "non-partonic" 
field. Strictly speaking, (two-particle) twist-4 operators containing two "non-partonic" 
fields also exist and have to be included. They can mix into three-particle and also four- 
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particle operators as in Fig. [T]i,e. However, contributions of such operators can always 
be dispensed off using equations of motion. Therefore, the results presented in our work 
are in fact sufficient for writing down arbitrary QCD evolution equations to the twist- 
four accuracy and e.g. calculation of the spectrum of anomalous dimensions of arbitrary 
twist-four operators. The kernels are written for the renormalization of coordinate-space 
light-ray operators [38J built of chiral fields and are manifestly SL(2) invariant. We 
believe that this form is most suitable in practical applications. The results can easily 
be recast in momentum space, in the form of evolution equations for generalized parton 
distributions [39] . The application of our formalism to DIS will be presented elsewhere. 

The presentation is organized as follows. Sect. 2 is introductory; we explain some 
basic ideas and the coordinate-space formalism. In Sect. 3 formal definitions are given, 
and we specify the conformal operator basis that is used throughout this work. Sect. 4 
contains a summary of the evolution kernels for quasipartonic operators from Ref. [32] 
which we rewrite in our language. The generalization of these results to the 2 — > 2 and 
2 — > 3 mixing kernels for non-quasipartonic operators is considered in Sect. 4 and Sect. 5, 
respectively, where we explain our method on simple examples. The complete results 
for the 2^3 mixing kernels are collected in Sect. 7. The final Sect. 8 is reserved for 
conclusions. 



In this work we will use the light-ray operator formalism in the spinor representation. 
The basic elements and the notation are explained in what follows. 

2.1 Light-ray operators 

We refer to non-local gauge-invariant operators with all the fields lying on a light-like 
line n 2 = as light-ray operators [ID] . The simplest example is 



2 Background 




(2.1) 



where q(x) is the quark field and [-21,-22] is the (light-like) Wilson line 




(2.2) 



Here and below we use "plus" for the projection on the light-cone direction 



a+ = n^a 




and also a shorthand notation 



Z\2 — Z\ — Z2 




UZi + uz 2 , 



u—l — u. 
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Because of the light-like separation, the operator in (12.11) contains additional (ultraviolet) 
divergences apart from usual field renormalization, which have to be renormalized. For 
practical purposes it is sufficient to define the renormalized light-ray operator [0]r as 
the generating function for the renormalized local operators 

[o] R ( Zl ,z 2 ) = J2H kHN - fcv \- qD + 7+ D+ (2 - 5) 

Note that all local operators appearing on the r.h.s. of (I2.5P have the same (geometric) 
twist t = 2, defined as t =dimension-spin |41j . 

The scale-dependence of the renormalized light-ray operator is governed by the re- 
normalization-group (RG) equation 

d J a s 

where H is the integral operator [38] 



( J^- + fa) j. + |^h) [0(z u z 2 )) r = , (2.6) 



0](z 1} z 2 ) = 2C F { jf ^ [20(zi, «a) - a0(*? 2 , z 2 ) - aO(«i, 



da ^° d/3 , 4) - ^0^, z 2 ) } . (2.7) 

Expansion of Eq. ( 12.7ft in powers of zi, z 2 generates the mixing matrix for local operators 
(12. 5p . Staying with the generating function (light-ray operator) offers, however, several 
technical advantages. One of them is that H and hence the RG equation in this form are 
manifestly covariant under the SL(2, R) transformations of the light-cone coordinates 
which correspond to the collinear subgroup of the conformal group [51] , see Sec. 2.3. 
Another advantage is that the RG equation (12. 6p . (12. 7p is completely general. It can be 
rewritten as the evolution equation for the generalized (quark) parton distribution [39] 
and reduces to DGLAP and ERBL evolution equations in the appropriate kinematic 
limits. For example, in the case of the total cross section of deep-inelastic scattering 
there is no momentum transfer between the initial and the final state, so that light-ray 
operators that differ by a total translation can be identified: (P\G(z\ + y,z 2 + y)\P) = 
(P\0(z h z 2 )\P). In this case Eq. (T22D simplifies to [38] 

[H • O] (0, z) — —2C p f duK(u)O(0,uz) (2.8) 

Jo 



with the kernel 
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Writing the DIS matrix element as a Fourier transform of the parton distribution 

(P\O(0,z)\P) = 2P + J dxe- lP+zx F{x) (2.10) 
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and picking up, e.g. the quark contribution, F q {x) = F(x),x > 0, Eq. (12.61) becomes 

f±. Fq {x^) = ^C F f ^-K(y)F q (x/y,n) (2.11) 
afi Z7r j x y 

so that CfK(x) is nothing but the familiar DGLAP quark splitting function P q ^ q . 
On the other hand, moments of K(x) correspond to anomalous dimensions of (flavor- 
nonsinglet) local operators: 

ln = -2C F [ dxx n ~ 1 K(x) = C F 
Jo 

The expressions in Eqs. (12. 6p . (12. 7p give an example of what we are aiming at for the 
description of higher-twist operators. The advantage of this, coordinate space formula- 
tion is that conformal symmetry is manifest and also there is an immediate relation both 
to renormalization of local operators and DGLAP-type equations for the corresponding 
(multi) parton distributions. 

2.2 Spinor Representation 

We use spinor formalism and follow the conventions adopted in Ref. [SB]. To this end, 
each co variant four- vector x^ is mapped to a hermitian matrix x: 

where <t m = (1, a), <r M = (1, —a) and a are the usual Pauli matrices. The Dirac (quark) 
spinor q is written as 

q=(x? > ,$«), (2.13) 

where ip a , x 13 are two-component Weyl spinors, ip a = (ip a ) , X a — (x°) ■ The gluon 
strength tensor F^ can be decomposed as 

F af3,af3 = a aa a pp F HV = 2 ( € a$fa(B ~ ^apf&fi) ' ( 2 ' 14 ) 

Here f a p and f & g are chiral and antichiral symmetric tensors, /* = /, which belong to 
(1, 0) and (0, 1) representations of the Lorenz group, respectively. For the dual strength 
tensor F^ v = \e^ vpa F pa one obtains 

iF af3,dP = 2 ( e a$f^ + taflfap) > ( 2 -!5) 

so that f al 3 (f & g) can also be identified as the selfdual (anti-selfdual) component of the 
field strength. The corresponding explicit expressions are 

fop = \ {D a a A aP + DfA aa ) , f. $ = i (D a a A a$ + D«A aa ) , (2.16) 



4(^(n + l)+ 7s ) 



n(n + 1) 



. (2.12) 
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where the covariant derivative is defined as = d^ — igA^. For convenience, we present 
the expressions for Dirac matrices in the spinor basis: 



7" = ' , c~ = ^ » , 7, = ^ " ^ I ■ (2.17) 

Here a My = §[7^,7"], 75 = i7°7 1 7 2 7 3 and 

(an J = \[o*o» - , {*n% = £[w - ■ (2.18) 

Any light-like vector can be represented as a product of two spinors. We introduce two 
independent light-like vectors as 

f^aa X a \a , Ti , 

n a a = Hafia, h 2 = , (2.19) 

where A = A*, /x = fj,<. The basis vectors in the plane transverse to n,n can be chosen 
as fi a Xa and X a p>a- An arbitrary four- vector be represented as 

X a a = Z X a \a + Z fi a p>a + W Xafia + W fi a Aq. , X 2 = (/iA) (A/x) — w] , (2.20) 

where z and 5 are real and w, w = w* complex coordinates in the two light-like directions 
and the transverse plane, respectively. 

The "+" and "— " fields are defined as the projections onto A and /x spinors, respec- 
tively: 



i> + = A a ^« , 


X+ 


=X a X a, 


/++ = 




i> + = A a ^d , 


X+ 




/++ = 




if)- = /x°V« , 




- = jF^a 


/+- = 




etc. Fields with free spinor indices 


can 


be written 


in terms of the 


"+" and 



(2.21) 



nents, e.g. 

(fl\)lp a (z) = \a4>-(z) - fl a 1p + (z) , 
(^X) 2 fa(j(z) = XaXp f—(z) - (\ a flp + A/J/X a ) f+-(z) + fl a fi/3 f++(z) ■ (2.22) 

Note that in difference to Ref. [36J we do not impose any particular normalization condi- 
tion on the auxiliary spinors A, \i. Without loss of generality one can put e.g. (/xA) = 1. 
However, keeping this factor proves to be convenient as it allows one to keep track of the 
balance of A, /x, A, /x in the equations. 

The light-ray operator (12.1 j) is decomposed as 

^1)7+^2) = ^+(zi)^+(z 2 ) + x+{zi)x+{z 2 ) , 
g(^i)757+g(^) = tMzi)tMz 2 ) - X+{zi)x+{z 2 ) , (2.23) 
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where q(z 2 ) = q(z 2 n) etc., and the Wilson lines are implied. 

As a less trivial example, consider the three-particle light-ray operators [38] 

Sp(zi,Z2,Za)=gq(z 1 ) F li+ (z 2 ) ±F li+ (z 2 )i^ l+q(z 3 ) (2.24) 

which contribute to polarized deep-inelastic scattering to twist-three accuracy (to the 
structure function g 2 (x, Q 2 ))- Going over to = cr^S^ one finds easily 

S+ & =2g [X^ + ( Zl )f +a (z 2 )^ + (z 3 ) + A a x+(2i)/+«(22)x+(^)] , (2-25) 

S aa = 2 9 [^aX+(zi)f+a(z 2 )x+(z 3 ) + K4>+ fcl ) f+& (z 2 )^J + (z 3 )] . (2.26) 

The operators contain both twist-three and twist-four contributions. Indeed, expan- 
sion of ip + (z\) f +a (z 2 )ip + (z 3 ) (and similarly for the other terms) at short distances goes 
over local operators O a ,ai,...a fc+1 ,d l ,...,d* (where k is the total number of covariant deriva- 
tives) which are symmetric in all dotted indices di, . . . , d& and also in the subset of the 
undotted ones, a±, . . . , aik+i- m order to separate the (leading) twist-three contribution 
one has to symmetrize in the remaining undotted index a. The answer can be written 
as follows (cf. [38]): 

[^i)/ +a (^ + (^3)] tw - 3 = ^J drr^^U^rz^Mrzs) • (2.27) 

As it should be, the twist-three part of this operator involves only "+" fields. For higher- 
twist operators twist separation becomes rather cumbersome in the "vector" formalism 
and going over to the spinor basis yields considerable advantages. We will encounter 
further examples in what follows. 

2.3 SX(2,R) invariance 

It is well known that conformal symmetry of the QCD Lagrangian imposes nontrivial 
constraints on the structure of RG equations to the one-loop accuracy, see e.g. Ref. [34|. 
For fields "living" on the light-cone, $(^n), it is sufficient to consider the so-called 
collinear subgroup SL(2, IR) of the conformal group 5*0(4, 2), corresponding to projective 
(Mobius) transformations of the line x = zn: 

z — > — ^— , ab — cd = 1 , 

cz + a 

where a, b, c, d are real numbers. A field with definite spin projection s on the light-cone 
transforms according to the irreducible representation of the SL(2, M) group with the 
conformal spin 

j = i(£ can + s) = £ can — E/2 , (2.28) 
2 



s 



where £ can is the (canonical) dimension and E = £ ca,n — s is the collinear twist. Action 
of the SL(2) generators on quantum fields can be traded for the differential operators 
acting on field coordinates. In this representation the generators become 

S + = z 2 d z + 2jz, S = zd z +j, S_ = -d z . (2.29) 

They obey the standard commutation relations 

{S + ,S_]=2S , [S ,S ± ] = ±S ± . (2.30) 

A finite form of the group transformations is 

Note that the field decomposition in "+" and "— " components as in Eq. f)2.22p is equiv- 
alent to the separation of different spin projections, s = +1/2 and s = —1/2, resulting 
in different values of the conformal spin j = 1 and j = 1/2 for ip+(z) and ip-(z), respec- 
tively. Similarly for gluon fields: /++ corresponds to j = 3/2, /+_ to j = 1, etc. 

The functional form of the "Hamiltonians" EI that appear in the RG equations for 
light-ray operators (12. 6p . cf. Eq. (12.71) . is constrained by the SL(2) invariance. For 
illustration, consider the simplest case: A light-ray operator built of fields with conformal 
spins ji and j 2 mixing into a light-ray operator built of the fields with the same spins. 
We are looking for an invariant kernel 7i\ 2 acting on functions of two variables (p(zi, Z2) 
which transform according to the representation T J1 <g> T- 72 : 

V ( Zl , z 2 ) - [I*(g) ® T»{g) V ]{ Zl , z 2 ) = {cz x + d)-^(cz 2 + d)~ 2 » V (^^, 

\czi + d cz 2 + d / 

(2.32) 

The SL(2) invariance means that H commutes with group transformations: 

Hv2 T n (g) ® T»(g) = T*(g) ® T j2 (g) H 12 . (2.33) 

There are different ways to define an operator. The first one is to specify the eigen- 
values of 7ii 2 on each irreducible component in the tensor product decomposition 

00 00 
T jl ® T J2 = © T^+ n , <p( Zl , z 2 ) = <Pn(zi, *2) • (2.34) 

n=0 n=0 

Let 

Tin <Pn{zi, z l) = Kifnizt, z 2 ) . (2.35) 

On the other hand 

S 2 2 Vn(z 1} z 2 ) = j n (j n - l)<p n (zi, z 2 ) , j n = n + j x + j 2 , (2.36) 
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where 



Sf 2 = -d x d 2 zf 2 + 2{ji - l)d 2 z 21 + 2(j 2 - l)d lZl2 + (ji + h - + i2 - 2) (2.37) 

is the two-particle Casimir operator, S 2 2 = (Si + S 2 ) 2 . Here Si and S 2 are the generators 
(I2.29P acting on z\ and Z2 coordinates, respectively, dk = d/dzk and Z12 = 21 — 22, cf. 
(12.41) . It proves to be convenient to define the operator J\ 2 as a formal solution of the 
operator equation 

£12 = Jn{J\2 - 1) , Jl2<fn{zi, Z 2 ) = j n <f n {Zl, Z 2 ) ■ (2.38) 

Since H12 and J% 2 have the same eigenf unctions, they can be diagonalized simultaneously 
which means that TC\ 2 = h{J\ 2 ) where h(j) is an ordinary function which can be found 
by expressing the eigenvalues h n (I2.35H in terms of j n (12.381) . For example, the operator 
in Eq. (12.71) can be written as 



H = 2C f 



^(J 12 + l) + ^(Ji2-l)-2^(l)-^ 



(2.39) 



where ip(x) = dlnT(x)/dx is the Euler ^-function. This is the most concise form. 
Moreover, as noticed in Ref . [37] , the renormalization group kernels involving fields with 
other spin projections can be obtained from this expression replacing the SL(2) Casimir 
operators by the ones of the full four- dimensional conformal group 0(4, 2). We will use 
this technique in Sect. 5. 

Often it is preferable to have a more functional definition of an operator. It can be 
shown that action of any SL(2, R) invariant operator 7ii 2 on a function (p(zi, z 2 ) can be 
written in the form 



[WiadCsi,**) = / daj^ dPa 2 K- 2 p 2 »- 2 u(^£) v^ 2 ,4), (2-40) 



where u(x) is an arbitrary function of one variable. For the same example in Eq. (12. 7p 

[H ■ 0}( Zl , z 2 ) = 2C F J\a £dp {5 [0( Zl , z 2 ) - 0(z? 2 , 4)] - 0(*? 2 , 4) " s}. 

(2.41) 

We stress that the conformal symmetry only becomes manifest in the properly chosen 
operator basis. Problem is that light-ray operators involving "— " components of the fields 
can mix in operators containing "— " or transverse derivatives that do not have, in general, 
good properties under conformal transformations. A general approach how to deal with 
such contributions in a conformally covariant way was developed in Ref. [36]. The 
idea is to expand the light-ray operator basis by adding primary fields with a particular 
transverse derivative, e.g. D^ + ip + , and at the same time eliminating the other derivative, 
D + -ip + , using equations of motion. We will elaborate on this proposal in the following 
Section. 
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3 Conformal Operator Basis 
3.1 Light-ray fields 

We define one-particle quark, antiquark and gluon light-ray operators (alias light-ray 
fields) including the Wilson line 



[0, jz]<&(z) = Pexpj-zgz J du(n ■ A)(itzn)j $(zn) . (3.1) 

The Wilson lines are always assumed although in many cases they will not be written 
explicitly. Here and below we use a shorthand notation $(z) for $(nz) etc. 

Light-ray fields can be viewed as generating functions for local operators with covari- 
ant derivatives that arise through the (formal) Taylor expansion 

[0, zMz) =]T^(n- D) k m = ^n(*D\) k m . (3.2) 

k k 

Note that all local operators on the r.h.s. of (13.21) have the same collinear twist as the 
field $ itself since each (n-D) derivative adds one unit of dimension and spin projection, 
simultaneously. 

It is implied that the Wilson line in Eq. fl3.ll) is written in the appropriate repre- 
sentation of the color group. In order to unify the notation we introduce the SU(N) 
generators acting on quark, antiquark and gluon fields as follows: 

= T^,ip 1 ' for the quark fields ip, x 
t a $ = { (t a ipy = -T^ft' for the antiquark fields x ( 3 - 3 ) 
_ ^bab'jb' f or ( a nti) self-dual gluon fields /, / 




where T a are the usual generators in the fundamental representation and i abc are the 
structure constants. Then, in particular 

Dip = (d — igA a T a )ip , 

Dip = ip(d +igA a T a ) . (3.4) 

We tacitly assume existence of rip quark flavors; flavor indices will not be shown explicitly 
in most cases. 

We will use the notation $ + for the chiral "+" fields and $ + for the antichiral ones: 

$+ = {iP + ,x+J ++ }, = {ip + ,X+J++}, 

= {^_, x-, /+-} , = X-, /+-} • (3-5) 

The fields $ + ,$ + have collinear twist E = 1, while have E = 2. In addition, 

one-particle E = 2 operators can be constructed by adding transverse derivatives to 
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/++ 




/+- 






3 


1 


3/2 


1/2 


1 


3/2 


2 


E 


1 


1 


2 


2 


2 


2 


H 


1/2 


1 


-1/2 





3/2 


2 



Table 1: The conformal spin j, collinear twist E and helicity H of the light-ray primary fields. 
The quantum numbers of the x+ an d X- fields coincide with those for tp + and ij}—, respectively. 
The anti-chiral fields -0+ , . . . have the same conformal spin as their chiral counterparts and the 
opposite chirality. 

E — 1 "plus" fields, e.g. D^ + ip + , D + ^ip + , where _D_ + = (/IDA) etc. Some of them 
(but not all, in a general situation) can be excluded from consideration with the help 
of equations of motion (EOM). A systematic procedure to treat such contributions is 
presented in Ref. where we have shown that EOM can be used to eliminate a "half" 
of the transverse derivatives in such a way that the remaining fields (with a derivative) 
transform as primary fields under the SL(2,M) transformations. In the above example 

one has to eliminate ip + in favor of D- + ip + . 

The complete basis of light-ray one-particle primary fields that is sufficient for the 
studies of twist-4 operators includes the following operators: 

X = {$+ , <J + , , (3.6) 

Each primary field X carries two more quantum numbers in addition to the conformal 
spin — collinear twist E and chirality H — which are eigenvalues of the two remain- 
ing generators of the full conformal group that commute with the light-cone SL(2, R) 
collinear subgroup, see Table. [TJ The fields t/>± and x± have the same quantum numbers, 
so that we display them for the ^-field only. The anti-chiral fields <£> have the same 
conformal spin as the chiral ones $ but opposite chirality H — > —H. 

3.2 Composite light-ray operators 

Gauge-invariant iV-particle light-ray operators can be defined as a product of primary 
fields 

0(z u ...,z N ) = SX( Zl ) <g> X(z 2 ) ® . . . ® X(z N ) 

= S h „ 4N ([0, z x \X{z y )f ([0, z 2 ]X{z 2 )) 12 . . . ([0, z N ]X{z N )Y N , (3.7) 

where X(z k ) are the fields from the set in Eq. (13. 6p . i\, . . . , are color indices and S , il ...j JV 
an invariant color tensor such that 

[(fx)JU + it 2 ) a k2l2 + ... + (t N ) a kNiN }S k> ... >iN = . (3.8) 

Here it is implied that the generators t a are taken in the appropriate representation, cf. 
Eq. (13.31) . The condition in Eq. (13. 8p ensures that 0(zi, . . . , z^) is a color singlet. 
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Composite light-ray operators (13. 7p transform as a product of primary fields under 
collinear conformal transformations 



k 



o{ ZX , . . . , z k ) - n( CZi + d)~^ o( Z [, . . . , 4) , 



(3.9) 



i=i 



where z ; = (az + b)/ (cz + <i) and ji is the conformal spin of the i-th "constituent" field. 

The collinear twist of a multiparticle light-ray operator is equal, obviously, to the 
sum of twists of the fields, E = E\ + . . . + En- For a given N, the lowest possible 
twist is E = N and it corresponds to operators built of "plus" components of the fields 
only, X = {$+, $+}• Such operators are known as quasipartonic [32]. The structure of 
logarithmic ultraviolet singularities in the relevant one-loop Feynman diagrams is such 
that the operators with less fields can mix into the operators containing more fields, but 
not the other way around. As a consequence, the set of quasipartonic operators is closed 
under renormalization. The corresponding RG equations were derived in Ref. [32]. They 
are sufficient for a calculation of the scale-dependence of arbitrary twist-three observables 
in QCD, e.g. in polarized DIS and also leading-twist baryon distribution amplitudes. 

Multiparticle operators with E = N + 1 built of one "minus" field with E = 2 and 
N — 1 "plus" fields with E — 1 are the next in complexity. They are subject of this 
paper. These operators can mix among themselves and also with N + 1-particle quasi- 
partonic operators. The RG equations for such operators derived below are sufficient for 
a calculation of the scale dependence of arbitrary twist-four observables. 

As an illustration, consider a typical twist-4 operator iqF^n^^j^q which contributes 
to the DIS structure functions at the level of power-suppressed l/Q 2 corrections in the 
"longitudinal" operator basis of Ref. [23]. Going over to spinor notation one obtains 



Each of the four terms can now be rewritten as a combination of two "plus" fields and 
one "minus" field, e.g. 



Twist-four operators built of two "minus" fields also exist, e.g. but they can 

(and should) be eliminated using EOM so that one does not need to consider them 
explicitly. In this particular case (DIS) twist-four operators involving Z)_ + $ + , D_ + <& + 
can be eliminated using EOM as well, however, only at the cost of loosing manifest 
SX(2,R) covariance of the evolution equations, see Ref. [36]. To avoid confusion, we 

stress that the operators with transverse derivatives D ^ + , _D_ + <i> + are introduced 

here in order to maintain conformal covariance as an extension of the "longitudinal" 
operator basis of Ref. [23] • They are not the same as the "transverse" operators of Ellis 
Furmanski-Petronzio (EFP) [23] that are advantageous in another aspect: they yield 
simpler coefficient functions. 



iqF^n^jsQ =^+(f^)+ ~ 0/)+^+ - X+(fx)+ + (xf)+X+ ■ 



(3.10) 



(3.11) 



13 



3.3 Renormalization-group equations for light-ray operators 

Operators with the same quantum numbers mix under renormalization. Let Oi{X) 
i — 1, . . . , L be the complete set of such operators. A renormalized operator is written 

as 

[O l (X)] R = Z tk O k (X ), (3.12) 
where X = Z X X is the bare field. Renormalized operators satisfy the RG equation 

(3.13) 

Here (3(g) is the (QCD) beta function and 



d ™ „ 1 
dfj, 

is the matrix of anomalous dimensions. To the one-loop accuracy in dimensional regu- 
larization D = 4 — 2e one obtains 



Z = 1 + 



4ne 



■H 



and 



ot s ™ 
7 = — H. 



(3.14) 



The operator HI (Hamiltonian) has a block-triangular form (at one loop). It follows from 
the fact that the N— particle operators can only mix with M > N— particle operators, 
thus 



0(3) 

) 



/ H (2-2) H (2- > 3) . . . 

o rf 3 ^ 3 ) rf 3 ^ 4 ) 



0(3) 



o u {N ^ ) 



(3.15) 



Further, it follows from the inspection of Feynman diagrams that the diagonal blocks 
are given by the sum of two-particles kernels, M ik 



(2-2) 
ik 



i,k 



The general structure of the kernels is 

M^lX^( Zl )®X-(z 2 )] = ^^^gtWgl 1 



1 i\ i 



'A 



X^](z u z 2 ). 



(3.16) 



(3.17) 



Here [C g ]*-^ 2 is a color tensor, is an SL(2,M) invariant operator which acts on 
coordinates of the fields, and q enumerates different structures. Except for the cases 
when X <g) X — ip <g) ip, X®Xi f ® /, the operator EI does not change the components 
of the primary fields, i.e. H. : tp <g> ip — > ip <8> ip,ip <8> f — > ip <8> f and so on. 
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Similarly, for MS n ~~* n+1 ^ one gets 



H («-n+i) = J^Wk ( 3 - 18 ) 

i.k 

with 

uf 2 ^\x^{ Zl )®X^{z 2 )] = EE^E!^* 1 ^ 1 ^^]^.^), (3-19) 



where the (2 — > 3) kernels Tt^ are of course different from the (2 — ► 2) ones in Eq. (I3.17p . 
and so on. 

The two-particle kernels H^ 2 ' ) , H^ - * 3 " 1 correspond to the counterterms to the product 
of light-ray fields (13.11) with open color indices M = X(zi) <S>X(z 2 ), where X(z) belongs 
to the set (13. 6ft . Such objects are, obviously, not gauge invariant and discussing their 
properties one has to specify the gauge fixing scheme. We will use the light-cone gauge 
and the background field formalism (see jl3] for a review). 

We recall that in this formalism one splits the fields in the quantum and classical 
components, q —>■ q d + q, — > A d + a M . Taking into account quantum corrections to 
the composite operator M = M(q, A) corresponds to the calculation of the path integral 
over the quantum fields 



((M(q d 1 A d ))) 



J VqVa M{q d + q, A d + a) exp S^S R (q d + q, A d + a) + ^a 2 + J . 



(3.20) 



Here Sn(q, A) is the renormalized QCD action, a + = (a ■ n), where n 2 = and £ is the 
gauge fixing parameter which, in the light-cone gauge, has to be sent to zero, £ — > 0. One 
has a certain freedom to choose the transformation properties of quantum and classical 
gauge fields under gauge transformations [13]. In particular one can assume that 

a' = UaU^ , A' = UAU^ - -Udrf . (3.21) 

9 

The light-ray fields ( 13.11) transform homogeneously, X' = U(z = 0)X, so that M(q,A) 
transforms as the product of gauge matrices, U(0) <S> U(0), taken at space-time point zero 
and in the appropriate representation of the color group. As it is easy to see, Eq. ( 13.211) 
guarantees that ((M(q cl , A d ))) transforms in exactly the same way. This means that 
the counterterms to the product X(z\) ® X(z 2 ) are given in terms of products of the 
fields X themselves. For example, the gluon field A d can only appear as a part of the 
field strength tensor, F d , or inside a covariant derivative. 

The advantage of using the light-cone gauge is that it makes explicit the two-particle 
structure of the RG equations. The price one pays for this property is, however, breaking 
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of the Lorentz invariance. In particular, "+" and "— " components of the fields are 
renormalized in a different way 



± o 



Z±q±, 



[A 



1*1(0) 



At one loop 



and 



with 



1 + ^C F 
8vre 



R 



A/2 



87TC 



-c f 



\nn) 



(3.22) 



(3.23) 



1 + 



a s 

47T6 



11 AT 2 

—N c n f 

3 3 ; 



I 



2ne' 



(3.24) 



The renormalized coupling constant is related to the bare one as go = Z^ l ^ 2 g. In physical 
quantities such as the S'-matrix or correlation functions of gauge-invariant operators, the 
Lorentz invariance is, of course, restored. 

For more details on the renormalization of QCD in the light-cone gauge see Ref. [SI 
and references therein. 



4 Quasipartonic Operators 

The primary field X (13. 61) has six different "plus" components so that one has to know 21 
(2 — > 2) kernels to describe the renormalization of an arbitrary quasipartonic operator. 
Parity and charge conjugation symmetry leave only seven independent kernels which 
were all calculated in Ref. |32j. These results are summarized below. 

Conformal symmetry dictates that the two-particle kernels acting on the coordinates 
of the fields must have the general form in Eq. (12.401) 0. In one-loop calculations in QCD 

* Group theory tells us that a nontrivial kernel Ti. which satisfies (|2.33j) exists if and only if j\ + j<x = 

j[ + j' 2 + m, where m is an integer number. For our purpose it is sufficient to consider the following 

cases: j\ — j[, ]2 — 3%, ji = j[ ± 1/2, ji = j' 2 ± 1/2. All of them can be reduced to the first one using 
that ^'1^^1-1/2^2-1/2 = •.; an d ^132^1+1/2,32+1/2 = z -i W ih^hh _ 
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only a limited amount of "standard" functions appear: 

" da 



[H<p]{z u z*)= / — 2 V (z 1 ,z 2 )-a 2j i- 1 p(z? 2 ,z 2 )-a 2j >- 1 <p(z 1 ,z% 1 ) , (4.1) 

J a I J 

[H d <p](z 1 ,z 2 ) = I daa^-W^^z^z^), (4.2) 
Jo 

[H + V ]{z u z 2 ) = f da r^a^- 2 ^-V(^ 2 >4)> (4-3) 
Jo Jo 

[n + ¥>](zi,z 2 ) = j\a£dPa^- 2 p^- 2 f^j ^ 2 ,4), (4-4) 



[H-y\{z u z 2 ) = / 1 da FdPaP- 2 ^- 2 (4.5) 



all of which correspond to a "diagonal" mapping of conformal spins T J1 ®T^ 2 — > T J1 ®T^ 2 . 
In particular the operator B(2ji, 2j 2 )~ 1 H d , where B(x, y) is Euler beta function, defines 
the projector P onto the invariant subspace T^ 1 ^ 2 with the lowest spin n = in the 
tensor product decomposition in Eq. (12.341) . We will also use the notation n for the 
projector 

n = / - v = i - B{2j u 2j 2 y 1 n d . (4.6) 

More functions, Ti, e '^ with < k < 2ji, are needed for the case that the conformal spins 
are reshuffled as T jl <g> T j2 -> T jl " k / 2 ® T j2+k / 2 

[H e 1 f V ](z 1 ,z 2 )= C daa 2 ^- k - x a k -\{z a X2 ,z 2 ), < k < 2j x . (4.7) 
Jo 

Last but not least, we need two constants 

3 

^ = 4. 

a 9 = 6 /4AT c , b = yiV c - -n/ , (4.8) 

which correspond to the "plus" quark field (13.221) and transverse gluon field renormal- 
ization (I3.24[) in the axial gauge 

a a 
Z q = l + ^a q C F , Z g = l + -^a g C A . (4.9) 

Thanks to the constraint in Eq. (13. 8p the contribution due to field renormalization for 
the color-singlet operators ( 13. 71) can be rewritten as [32] 

v 1 - ^ E a x* C x) S ^ ® ■ ■ ■ ® = 

i 

= S ( 1 + 2^ E(°* + ® *i) ( X i ® • • • ® X ^v) , (4.10) 
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where cr^ = {cqj&g} an d Cx* = {Cf,Ca} for quarks and gluons, respectively. This 
property allows one to include the field renormalization factors in the definition of the 
two-particle kernels [32] . 



4.1 Coordinate-space representation 



A) X ® X = <g> if>+, ^+ ® x+, ^+ ® ^+ ® x+, X+ ® X+, X+ ® X+} 



In compact notation that we will also use for the other cases 



mx( Zl )®x(z2) 



-2(t a ®t a ) 



H-2o n 



X{zx)®X{z 2 ) 



(4.11) 



The operator in square brackets on the r.h.s. acts on the function (p(z\, z 2 ) = 
X(zi) <S> X(z 2 ) according to Eq. (14.51) where one has to substitute the values ji = 
j 2 = 1 for the conformal spins of the participating fields. Here and below a symbolic 
expression X <g> X implies that both fields have open color indices and 

{t b ® t b ) X{z x ) ® X(z 2 ) = t b X{ Zl ) ® t b X{z 2 ) = t^X^zt) t b i2i ,X^{z 2 ) . 



B) X <g> X = {</>+ ® x+, 4>+ ® x+> ^+ ® X+ ® X+} 



For quarks of different flavor: 

H X(zi) ® X(z 2 ) = -2(t b ®t h ) 



n-n 



2a n 



X(z 1 )®X{z 2 ). 



(4.12) 



The generators t a acting on the quark and the antiquark fields are defined in 
Eq. (13. 3p . After the contraction of open color indices S ili2 = S ili2 , cf. (13 .7p . this 
expression reproduces the result in Eq. (12.71) . 

For the quark- antiquark pair of the same flavor, X ® X = ip <S> ^, x ® X, there are 
two extra terms: 



wx\ Zl ) ® x\z 2 ) = ...- 44 n d r{z x , z 2 ) 



-2iz l2 [{t a t% 



+ 2(t b t 



flJzJ^Rjz,). (4.13) 



n + + n 

The ellipses stand for the contribution of Eq. (I4.12p . Here and below 

Pv2P(zi,z 2 ) = (p(z 2 ,z 1 ) (4.14) 
is the permutation operator acting on the field coordinates and 

J a (z u z 2 ) = (^(zi)T a ^(z 2 ) + xi(zi)T a xi(z 2 )) , (4.15) 



where the sum runs over all possible flavors. The generators t a in Eq. (I4.13P have 
to be taken in the quark, = T^-, and the antiquark, = —T^, representations 
for the pairs ip l ® ^ and x % ® X 1 '■> respectively. 
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c) x a ®x = {f a ++ ®i, +1 ft + ®X+} 



UX a { Zl )®X\z 2 ) = -2{t b aa ,®t\ v ) 



H — (J q — (Jg 



2(t a 't% i ,P 12 H e d 1) X a '(z 1 )®X i '(z 2 ). 



X a \z x ) ®X i \z 2 ) 



(4.16) 



Here j\ = 3/2, j 2 = 1 and the SU(N) generator acting on the first (gluon) field 
(in the first line) is in the adjoint representation. 



D) X a <g> X = {/« ® V>+, /« <g> x+, fl+ ® V+, /++ ® X+} 



— 27i — a q — a g 



mX a {z 1 )®X\z 2 ) = -2{t h aal ®t b ll , 

+ A(t a 't a ) ii m- x a ' ( Zl ) ® x r (z 2 ) 



X a '(z 1 )®X i '(z 2 ) 

(4.17) 



E) X ® X = {/ ++ ® /++, / ++ ® / ++ } 



H X(^) ® X{z 2 ) = -2{t b ® t b ) 



H-2a R 



X{z x )®X{z 2 ). (4.18) 



F) X®X = {f ++ ®J ++ } 



MfZ + ( Zl )®f c ++ (z 2 ) = -2(t b aa ,®t b cc MH-M + -2H+ -2a g /^(zO ® £' + (z 2 ) 



+ l2{t a 't a ) cc ,H-fi + { Zl ) ® /;' + (^ 2 ) + — | 2ft + P 12 - P ac ) ITo J ac (z 2 , z x ) , (4.19) 

Zl2 



where 



J-(z 1 ^ 2 ) = ^^(^i)T a T c ^(z 2 )- X ^(^)T c T^(^i) (4.20) 



and 



The operator IIq 



PacJ aC (z 1 ,Z 2 ) = J Ca (z 1 ,Z 2 ) 



(4.21) 



n <p( Zl , z 2 ) = (I — 6WX^, z 2 ) = 6 [ da aa [<p( Zl , z 2 ) - y{z« 2 , *° )] (4.22) 

Jo 

is the projector onto the subspace J\ 2 > 2 in the tensor product decomposi- 
tion f)2.34p . n = 0( J \2 > 2): The function i/j(zi, z 2 ) = Ho(p(zi, z 2 ) vanishes on the 
line z 2 = zi, ip(z, z) = 0. 
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4.2 Invariant representation 



As discussed in Sect. 2.3, calculating the eigenvalues of the two-particle kernels on the 
set of functions that correspond to irreducible representations of the SL(2, R) group one 
can write the kernels as functions of the two-particle conformal spin operator, J( J — 1) = 
(Si + S2) 2 , cf. Eq. (I2.38p . This representation will be useful to restore the two-particle 
kernels for non-quasipartonic operators. 



A) A <g> A = {ip + ® v+, ® X+, ^+ ® V>+ ® X+, X+ ® X+, X+ ® X+} 



H X[z x ) <8> X{z 2 ) = -4(t a ® t a ) V( J) - VK 1 ) - ^ ® • (4.23) 



B) A ® X = {</>+ <g> x+, 4>+ ® X+, ^+ ® X+ ® X+} 



For quarks of different flavor: 



HI(zi)®I(z 2 ) = 

- 2{t a ® t a ) [V>( J + 1) + ^( J - 1) - 2^(1) - 2aJ A(zx) <8> A(z 2 ) . (4.24) 
For quarks of the same flavor (ip ® ip ox x ® X om y) there are extra terms 



H X{z x ) ® X(z 2 ) = ...--*« 5 J>2 J a ( Zl , z 2 ) 



a+b \ 



2iz 12 (t a t 



2 
3 
1 



' y (J+l)(J-2) 



2(-l) J P. 
J(J 



) —^) /?+&) ® fXM) ■ ( 4 -25) 



c) x a ®x = {fi + ®^+, 7^+®^+, 



H A a (^) <g> X\z 2 ) = -21 (t b aa , ®t*,) ^(J + l/2) + t/>(J — 1/2) - 2^(1) - <r, - a 6 



+ (t a 't a 



' iy-1/2 
J -1/2 



>X a (zi) ® A* (z 2 ) 



(4.26) 



D) A a <g> A = {/• ® /» <g> x+, /i+ ® /:+ ® X+} 



H A a (z!)®A l (^) = -2<( (4'®4) ^ + 3/2)+V>(J-3/2)-2^(l)-<7,-o- fl 



+ (^)«y2(-i) J ^ r(j + 3/2) ^"'fr)®*^)- ( 4 - 27 ) 
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E) X ® X = {f ++ ® / ++ , /« + ® fl + } 



M X( Zl ) X{z 2 ) = -4(t fe <8> t fe ) V( J) - " ff g ® ^(22) 



(4.28) 



F) X®X = / ++ ®/ ++ 



H r + (^x) ® /?,(^ 2 ) = -2<^ (C ® 40 + 2) + ip{J - 2) - 2^(1) - 2a 



+ (t a 't a ) c M-i) jT{J 2) 



2z 

+ — 



r(j + 2) 

P a jQ{J>2)J ac {z 2 ,z 1 ). (4.29) 



Z12 V^(^-l) 



Note that with the exception of quark-gluon transitions (that are only possible for 
Xi ® X 2 = f ® / or flavor- singlet ip®ip.,x®X pairs) , the pair- wise Hamiltonians contain 
two color structures only, t c ® t c and (t a t a )ar. Moreover, apart from the contributions 
due to field renormalization, the corresponding SL(2, M) invariant kernels can be written 
in a universal form, as a function of one parameter: the difference of field helicities 



h — I hx 1 — hx 2 1 • 

The corresponding expressions are 

Hi = —2[ip(J + h)+ tj){J -h)- 2^(1) - a Xl - a x , 2 

n 2 = -2(-iy 



, h T(2h)T(J-h) 



r{j + h) 



(4.30) 



for the first, t c ® t c , and the second, (t a 't a )u>, color structures, respectively. 

It is interesting to note that only these two Hamiltonians appear in supersymmetric 
extensions of QCD, see Refs. [4"6| I4T]. 



5 Non-Quasipartonic Operators 

As noticed in Ref . [37] , knowledge of the two-particle evolution kernels for quasipartonic 
operators allows one to restore the similar kernels for non-quasipartonic operators by 
"submerging" the SX(2,R) subgroup into the full conformal group £0(4, 2). 

To understand the idea, recall that the power of conformal symmetry in the descrip- 
tion of renormalization of the quasipartonic operators is due to the fact that the tensor 
product of any two one-particle representations of the SL(2, R) group is decomposed into 
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irreducible components without multiplicities, cf. Eq. (12.34)) . Different irreducible com- 
ponents can be labeled by the eigenvalues of the single quadratic two-particle Casimir 
operator S\ 2 ( 12.381) . It follows that the evolution kernels EI and the Casimir operator 
S\ 2 have the same eigenfunctions and can be diagonalized simultaneously; hence H can 
be written as a function of Sf 2 , as exemplified in Eq. (I2.39p . This is what we call an 
invariant representation. It is easy to see that the non-degeneracy is of principle impor- 
tance: if two operators exist with the same transformation properties, they can mix in 
an arbitrary way. 

The full conformal group 50(4, 2) is more complicated and the same property does 
not hold for the tensor product of arbitrary representations: in general, there are non- 
trivial multiplicities. The observation of Ref. [37J is that the degeneracy does not occur, 
however, for the products of representations of the special type that we need for the 
construction of the QCD operators. 

It can be shown that the tensor product of two primary fields (13.61) is decomposed 
into irreducible components of SO(4, 2) without multiplicities. Schematically 



and each term (component) in the sum can be characterized uniquely by the value 
of the quadratic (two-particle) Casimir operator Cf 2 of the full conformal group (cf. 
Appendix A). Following the same line of argument one concludes that the evolution 
kernel (Hamiltonian) HI acting on an arbitrary two-particle light-ray operator must be a 
certain function of Cf 2 alone, which means that any operator from a particular irreducible 
component T^ n ^ is an eigenvector of EI with the same eigenvalue E n . On the other 
hand, each T^ n ^ contains a SX(2,R) invariant subspace corresponding to quasipartonic 
operators. Thus the eigenvalues E n on the subspaces Tr*™) , and therefore the operator H 
itself, can be restored from its dependence on Sf 2 in the quasipartonic sector. 

Using this strategy, we derive below a complete set of 2 — > 2 evolution kernels for 
arbitrary operators containing one "plus" and one "minus" field in the coordinate-space 
representation. 

5.1 Fields of the same chirality 

In this section we calculate the two-particle evolution kernels for the tensor product of 
two light-ray fields of the same chirality. 

With the restriction to collinear twist E = 3, there are four such operators: 



where &± are defined in Eq. (13.51) . The first two, $_(^i) ® Q + (z 2 ) and § + (zi) ® §_(z 2 ), 
have the same quantum numbers so they can mix under renormalization and the same 




(5.1) 



n 



® $+(z 2 ), ® ®-(z 2 ), 

$+Czi) ® D_ + ^+(z 2 ), D_ + <$> + ( Zl ) ® $ + (z 2 ) 



(5.2) 
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is true for the the second pair, <8> D_ + § + (z 2 ) and D- + § + (zi) ® $ + (z 2 ). On 

the contrary, light-ray operators with and without an extra transverse derivative, e.g. 
&-(zi)(B)Q + (z 2 ) and $ + (zx)<g)D- + & + (z 2 ), cannot mix because they have different helicity, 
cf. Tab. 1. As a consequence, the evolution kernel has a 2 x 2 matrix form, e.g. 



H 2 i H 22 



(5.3) 



and similar for the other pair. 

From the group theory point of view, the diagonal kernels H n and EI 22 are the 
SL(2, K) invariant operators which map T- 71 (g>T J ' 2 — * T jl ®T^ 2 . The off-diagonal kernels 
intertwine the representations T J2 ® T jl and T jl ® T- 72 , 



H12 : T J2 <g> T n 



H 21 : T J1 ® T J2 -»■ T J2 ® T ix 



Our task is to find explicit expressions for these kernels. 

In what follows we explain the method in detail on the particular example of the 
operator ip_ The other cases are similar so we only present the results. 

As the first step, one has to compare the action of 5*0(4,2) and SL(2,M) Casimir 
operators on the tensor product of two "plus" fields. Using Eqs. (lA.lOj) and (lA.llj) from 
Appendix |A] one finds 



C 2 



12 



s 2 
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ip+{zi) (g> Tp+(z 2 ) = -d x d 2 z\ 2 ip+(zi) ® ^+(^2) , (5.4) 



where we used the explicit expression for Sf 2 (12.371) (for the case j\ = j 2 = 1) to arrive 
at the last equality. Note that the action of SO (A, 2) and SL(2, M.) Casimir operators on 
the "plus" fields do not coincide, but differ by a constant (lA.lip . 

Thanks to the SL(2, R) invariance it is sufficient to consider the action of the Casimir 
operators on functions of two variables that are annihilated by the "step-down" operator 
(Si + <S£")</?(zi, z 2 ) = 0. This condition is nothing but shift-invar iance: <p(z%,z 2 ) = 



(p(zi — z 2 ) = (f(zi 2 ). It is easy to see that (fi n (zi,z 2 ) 
corresponding to the eigenvalue (n + 2)(n + 1): 



"12 



are the eigenfunctions 



did 2 z 



z \ 2 — 3n(jn 1) z 



12 ) 



j n = n + 2 . 



(5.5) 



The second step is to find the explicit expression for the Casimir operator acting on 
the ijj- ® ip + fields. It has a matrix form 



r 2 
^12 



(5.6) 



The entries a, 6, c, d can be found using Eq. (1A.9j) and explicit expressions for the con- 
formal generators fl A .ID . This is easy for the diagonal entries a and d, because the terms 
in the last two lines in Eq. flA.91) do not contribute. One obtains 



; 12 



:l/2j 2 =l) 



+ 



d 



O-i 



'12 



:lj 2 =l/2) 



(5.7) 
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The expression for S\ 2 is again given by Eq. (I2.37p . but with different values of the con- 
formal spins as compared to the "plus-plus" case. Hereafter we will not show the values 
of Ji,j2 explicitly, unless this can lead to a confusion. It is always implied that param- 
eters that enter the definitions of the operators are determined by quantum numbers of 
the fields they act on. 

The direct calculation of the off-diagonal elements b, c is rather cumbersome because 
in this case all terms in Eq. flA.9j) have to be taken into account. These entries can easily 
be found, however, with the help of the following trick. 

Let us apply the operation A(<9/<9//) to the both sides of (15.61) . Taking into account 
that by definition ip_ = fi a ip a one obtains 



'12 + 



1p+ ® "0H 



a b\ fij) + <g> ip + 
c d) 1-0+ ® ^+ 



The l.h.s. of this equality is given by Eq. (15.41) so that obviously 



b = S 



2,CJ'1=J2=1) 
12 



a = d 2 z 21 , 



S 



2,(ji=i2=i) 

12 



- d = dx Z12 ■ 



(5.8) 



(5.9) 



It is easy to check that the operators b and c intertwine the representations T 1 / 2 ® T 1 
and T 1 <g> T 1//2 , 



b : T 1 ®T 1/2 -5i> T 1/2 g r 1 
Collecting everything, we obtain 



C?2 + ^ 



S 2 12 + 1/4 d 2 z 21 



dxzu S 2 12 + 1/4 
Remarkably enough, one can represent this operator in the form 

3~ 



C 2 2 4 



I J^2 — 1 



where 



12 



^2-^21 

diz 12 



(5.10) 



(5.11) 



(5.12) 



The eigenfunctions of the operator J 12 have the form 



^(z 1 ,z 2 ) = z 



12 



1 
±1 



so that 



C 2 + 



^=Jn(jn -l)^n, j+ = 71 + 2 , j~ = (n + 1) . 



(5.13) 



(5.14) 
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As stated above [37] , the Hamiltonian EI must be a function of the full Casimir oper- 
ator of the conformal group. This means that they share the same set of eigenfunctions 
and the eigenvalues E n of EI are functions of j n . Moreover, the expressions for E n = E(j n ) 
must be the same for the both collinear twist E = 2 and E = 3 sectors: 



M- 



E=2 n 

n2 



E{jn)zl 



12 5 



E=3j,± 



n i 



(5.15) 



so that the function E(j) is already known: It is determined by the Hamiltonian for the 
quasipartonic operator ip + <g> see Eq. (I4.23j) : 



E(j) = 2[m - ^(1) - a q ] . 



(5.16) 



(We omitted here the color structure.) 

The remaining steps are purely technical. First, using (I5.15p . we can restore the 
action of EI on arbitrary shift-invariant polynomials. It is easy to see that the dimension 



of the polynomial is conserved, so we only need to consider functions of the form 
Using the standard decomposition EI = E + \ip+){ip+\ + E^\ip~){ip~\ one gets 



'12- 



'12 



hxx(n) hn(n)\ ( a 
h 2 i(n) h 22 (n)J \/3 



'12 5 



hi 



ik 



(5.17) 



so that 

hn(n) =h 22 (n) 



1 



hl2 (n) =h 21 (n) = - (E(j+) - E(j')) = ? _ { 



+ E{j~)) = *l>(n + 2) + V(n + 1) - 2V>(1) - 2a q , 
1 



(5.18) 



Once the action of EI on shift-invariant polynomials (which correspond to highest weights 
of the corresponding representations) is established, one can easily restore the explicit 
form of the SL(2, R) invariant operators in Eq. f !5.3j) from the requirement that 



in Mia 
EI21 EI22 



'12 



h n (n) h 12 (n] 
h 2 i{n) h 22 (n] 



'12 



(5.19) 



One obtains 



12 




■0+ ® 0- 



(5.20) 
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and 

'(£i£)-^'*(*rt'*%)(tity <5 - 2i) 

which is our final result. 

For all other cases the calculation is similar so that we only present the resulting 
expressions for the 5*0(4,2) Casimir operator and the Hamiltonians. Let 

&- ® r v Mih,„ „ ^ _ f®+ ® 



where ji,j2 and z x = j\ + 1/2, z 2 = j 2 — 1/2 are the conformal spin of the "upper" and the 
"lower" pair of fields. Below we assume that j\ < j 2 so that in the case of a quark-gluon 
pair the first field is always the quark and the second one is the gluon. The Casimir 
operator takes the form 

C? 2 + k X2 = J 12 ( J 12 - l) , (5.23) 



where K12 is defined in Appendix |Aj Eq. ( 1A.11I) and has to be calculated for the corre- 
sponding pair of "plus" fields $ + cg> $ + . The operator J 12 is defined as 

? C ' \T[ -1- f ^2i<9 2 + 2i 2 \ , * 

Here and below the upper sign corresponds to X + and the lower one to X_. 
The Hamiltonians for the "doublets" 



X 4 



V-®^ AA-®x+\ fx-®^ /7+-®/+ 
M) + ® tp-J ' ® x-y ' Vx+ ® V>-y ' u+ ® x-y ' u++ ® /+ 



A < ('£ . (t + ® D ~ +X+ ) , • • • , ® >, (5.25) 



® v^+y ' v°-+^+ ® x+y ' " ' ' \D-+u+ ® /++; 

(which can be viewed as the higher twist £7 = 3 descendants of the quasipartonic E = 2 



operators of type _A_ and _E_ ) take the form 



,e,l 



ELY ± (^, * 2 ) = -2(f ® t b ) [ H °?f 2 ~ ±Hl2 ) X ± ( Zl , z 2 ) . (5.26) 

±/t 21 ri-a Xl x 2/ 

Here 0x1^2 = +<^x 2 , &x = {c g , cr g } and we suppress all color indices (cf. Eq. (14. lip ). 
The Hamiltonian in the antichiral sector is exactly the same. 
Next, for the operator pairs 
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and 



X i _ a (z 1 ,z 2 ) 



ij\ ® [D„ + f ++ ) a \ fx% ® [D-+f ++ ] a 



(5.2? 



which are the descendants of the quasipartonic operators of type C , one gets 



,e,l 



MX™(z 1 ,z 2 )=-2{(ti,®t b aal ) 



T (t a 't a ) lll P 1 



H 

7~L — o Xl x 2 



12 



rL 21 ~r rL 21 

TH e £ 



Xi a '( Zl ,z 2 ). 



(5.29) 



We remind that all kernels depend on the conformal spins of the fields they are acting 
on. 



5.2 Fields of opposite chirality 

The two-particle evolution kernels for the tensor product of two light-ray fields of oppo- 
site chirality can be calculated along the similar lines. A difference is that e.g. $_(zi) ® 
$+(,22) cannot mix under renormalization with $ + (zi) ® $_(^) (since chirality is con- 
served) and the same is true for the pair $+(zi) ® -D_+$+(^ 2 ) and D_ + $ + (zi) <g>$ + (z 2 ). 
At the same time, mixing of the light-ray operators with and without an extra transverse 
derivative, e.g. $_(zi) ® ®+(z 2 ) and ® Z^_ + <l ) + (-22), is allowed. 

For this reason, instead of (I5.22p . we have to consider the "doublets" 

X ^ {zu Z2) = ($5 1 |[d ' (5 ' 30) 

where we take $ = {V>, x, /}, * = {V>, X, /}; ji, J2 and i x = j x + 1/2, % 2 = j 2 + 1/2 are 
the corresponding conformal spins. 

The 5*0(4, 2) Casimir operator acting on the the light-ray fields (15.301) can be repre- 
sented in the form (I5.23P with 

^=-01+^-1/2)1-^^), (5.31) 

where 

F 12 = -did 2 z 2X + 2(h - l)d 2 - 2{i 2 - l)<9i . (5.32) 
The eigenfunctions of J are 

^i,^)^- 1 ^ 1 ), n>0, (5.33) 
^-(z l7 z 2 )=z^ 1 h 1 ). ;/:>(). (5.3 L) 
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where a n = —n — 2(ji + j 2 ) + 1 and 

Ji/>± = T(n + h + h ~ 1/2)^ • (5.35) 
Thus, for the Casimir operator (I5.23j) one obtains 

[C ( i 2) + K 12 ) ^ = 3 ±Un - 1) ^ , in = n + h + h ± 1/2 , (5.36) 



where, as above, the value of K\ 2 (lA.llj) has to be calculated using the helicities of the 



corresponding "plus" fields $ + ® $ + , cf. Tab. 1. 

Using these expressions and proceeding along the same lines as in Sect. 5.1 we can 
restore the Hamiltonian in the coordinate space representation. 

First, consider the "doublets" that arise as descendants of the quasipartonic operators 



of type |_B 
X = 



i>- ® V+ _ A / x- ® V>+ _ \ ( ^-®x+ \ ( x-®x+ 

^+ ® \D„ + i> + ) ' \ X + ® \D_ + 4, + ) ' \^ + ® \D_ +X+ ) ' \X+ ® 

(5.37) 

We assume here that the quarks ^_ Cg> V>+ are in the flavor-nonsinglet state, so that there 
is no mixing with gluons. In this case one finds 

fU + n d -2a q z 2 iHt 2 
MX( Zl ,z 2 ) = -2(t b ®t b )l J_ u fi-2>H + -2a l A " ( ~^- ( '" ) ' 3$) 

where n is defined by Eq. (14. 6p . 

Descendants of quasipartonic operators of type _D_ fall in two classes that have to be 
considered separately. 

First, consider the operators 

The corresponding evolution kernel takes the form 

UX ia ( Zl ,z 2 ) = - 2{(4 <g> + (t a 't a )^H 2 }x iV (^,^ 2 ) . (5.40) 

The 2x2 matrix kernels Hi and H2 are given by the following expressions: 

l + H d -a qg z 21 Hf 2 \ ( -H d z 21 H u 

— n H - 37^+2 - o qg J ' 1 Pi 2 H 2 iU -3Hi 2 

(5.41) 

where a qg = a q + o g and n is defined by Eq. ( 14. 6 p . 
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The second set of operators is 

In this case one obtains 

MX ai ( Zl ,z 2 ) = - 2{(t b aa ,®$ i ,)U 1 + (t a 't a ) ii M 2 ^X a ' i '(z 1 ,z 2 ) (5.43) 



with 



Hi = 



'H-H+ + 2H d - a qg Z2i (nt 2 + Wf 2 ) 

^n W - 4W+ - 2W+ - <r Wy 



H 2 =( (5.44) 



Next, we consider the renormalization of // operators which are the descendants of 
quasipartonic operators of type F . We define 



/+- ® fl 



u — I — -in fb 



J [Z U Z2) - [i D _ + ^ +{zi)T a TbMz2) _ x+ ( Z2)7 *T>iD- + x + (z 1 )) ' [bAb) 
where in the second operator doublet the sum over flavors is implied. One obtains 
MG ab (z l} z 2 ) = - 2[{t c aa , ® t c w ) M x + (t c ba , ® t c ah ,) H 2 )G a ' b \z 1 , z 2 ) 

+ —hm 3 - pAUoJ^iz^z^ . (5.46) 
Here P a b is the operator of permutations in color space. The kernels Hj take the form 



Hi 



'ft - + m d - 2a 9 z 2 i (w+ + 2Ht 2 ) 

-i-n W - 6(W+ + W+) - 2a, 



^i^r 2 \ / -^2iWi 2 

^21 / \^21 

The projection operator n is defined as 

n = i - P , 
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where the operator P is the projector to the null subspace of the Casimir operator, 
C 2 12 F = 0. It has the form 



Po 



^2H d + \v l -3z 21 H d \ 



1, 
4' 



(5.48) 



21 



/ 



where V\ is the SL(2,R) projector to the subspace J = 5/2 on the tensor product 

T 1/2 0T 1 



[V 1 (p]{z 1 ,z 2 ) = iz 21 [ daa 2 a[(d 2 -2d l )<p}{z? 2 ,z? 2 ) 
Jo 



(5.49) 



Finally, we present the result for the quark-antiquark pair in the flavor-singlet state, case 



B . Let 



x AB ( Zl ,z 2 )- )t W x-®x B + 



'> A ® i> B 

\A x> l n „t;b / ' V /o, i 



^1,22) = i 



al 



(5.50) 



(5.51) 



where in the first operator A, B are flavor indices of the quark fields and in the second 
operator the sum over flavors is implied. One obtains 



WX AB {z 1} z 2 ) = -2(t a <g> t a ) Hi X AB (^, z 2 ) - 6 AB ^ P J"( : , ! 



- 5 AB 2iz 12 (t a t% [H 2 + P a6 H 3 ] G^, * a ) 
where G afe (2;i,2; 2 ) is defined in Eq. ( I5.45P and 



(5.52) 



Hi 



n + n d - 2a q 

z 2 i 



z 21 H 



12 



H - 2Hf 2 - 2(7, 



Ho 



-z 21 W 



12 



— w 12 n 7^ 2 + 37Y 12 

^21 



+ ' 



-z 21 H 



12 



— ^i 2 n 47Y 12 

.221 



(5.53) 



6 Mixing with Three-Particle Operators 

As it was already mentioned in Sect. 13.31 light-ray operators with a different number 
of constituents can mix with each other. To the leading order in strong coupling this 
mixing has a triangular form ( 13. 15ft so that at the level of kernels we need the 2-^3 
contributions corresponding to the mixing of one "plus" and one "minus" primary fields 
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with the three-particle quasipartonic operators. It seems at first sight that the diagonal 
and off-diagonal blocks in the matrix (13 . 1 5[) are independent and a separate calculation 
is necessary to fix the missing off-diagonal terms. Surprisingly enough this is not the 
case. We will show below that all 2 — > 3 kernels are completely determined by the 2-^2 
kernels in the quasipartonic sector. 

In short, the main idea is that the necessary relations are imposed by Lorentz sym- 
metry. For example, applying the generator of translations in transverse direction to the 
light cone to the "plus" quark field one obtains 

i[P^, = = 2«9+^_ + igA^+ + EOM , (6.1) 

so that the similar transformation applied to the renormalized quasipartonic light-ray 
operator yields 

® ip+(z 2 )] R = [d^ + {z x ) <g> ^+{z 2 )} R + foMzi) ® d^ + {z 2 )] R 
= 2d zl [4j^{z 1 )®^ + {z 2 )] R + 2d Z2 [t{j + {z 1 )®^{z 2 )] R 

+ wlApxip+izi) ® ip+(z 2 )] R + ig[ip + (zi) ® A^ x ip + (z 2 )] R 
+ EOM. (6.2) 

Eq. (16.21) is an operator identity which must be satisfied by the corresponding opera- 
tor counterterms order by order in perturbation theory. The operator on the l.h.s. is 
quasipartonic and its renormalization is governed by the corresponding BFLK kernel. 
As we will see below, application of the transverse derivative to the two-particle coun- 
terterm will generate a sum of contributions of two-particle and three-particle operators. 
The operators ip-(zi) <8> ip+(z 2 ) and ip+(zi) g) ip-(z 2 ) in the second line in Eq. (16.21) con- 
tain both two-particle counterterms that we have calculated in the previous Section and 
the three-particle counterterms that we do not know so far. Finally, the three-particle 
quark-quark-gluon operators in the third line in Eq. (16.21) are again quasipartonic and 
their renormalization is described in terms of the two-particle BFLK kernels. Thus 
Eq. (16. 2p provides one with a relation for the 2 — > 3 mixing kernels for the ip- <S> ip+ 
light-ray operators and the BFLK kernels. The question is whether this constraint is 
sufficient to determine the 2^3 kernels in principle, and how to solve it in practice. 
We have found two possibilities to proceed. 

The first one is to derive another constraint, applying Lorentz rotations M MM instead 
of P^x to the quasipartonic ip + operator. We are able to prove that, taken together, 
these two constraints determine the 2 — > 3 kernels uniquely. An obvious advantage of 
this approach is that it only relies on exact Lorentz symmetry and so it may also be 
applicable beyond one loop. The main disadvantage is that derivation of the constraint 
imposed by the M MAt transformation is more complicated because it also affects the axial 
gauge fixing term in the QCD Lagrangian. The corresponding contributions must be 
taken into account but, at least to the one- loop accuracy, are relatively simple because 
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of cancellations between different Feynman diagrams that have a structure typical for a 
Ward identity. 

The second possibility which we eventually found to be the most effective, is to 
study properties of the constraint equation corresponding to Eq. (16.2p under the collinear 
conformal transformations. As it stands, this equation is not SL(2) invariant, but, as 
we will explain, it can be separated in two SX(2)-invariant equations that are already 
sufficient to determine the 2^3 kernels of interest. To the one loop accuracy, both 
approaches are equivalent and produce identical results. Some of the results presented 
below have also been checked by the direct calculation of relevant Feynman diagrams. 

It is worthwhile to note that the same technique can be used to determine the 2-^2 
kernels as well. In this case, however, the approach based on the construction of the 
5*0(4,2) Casimir operator [37], Sect. 5, is clearly advantageous. 

In what follows we explain the details of our approach on the concrete example of 
the ip_ (g> tp + operator and then present the results for all kernels in question. For this 
example, we consider the following set of operators: 

0{ k (z 1 ,z 2 )=i/; i _(z 1 )®i/; k + (z 2 ), 

Of\z u z 2 , z 3 ) = ® ^2) ® f h ++ {z,) , (6.3) 

where i,k,b are the color indices that in most cases will not be displayed explicitly. 
The first operator in this list is quasipartonic. To one-loop accuracy the renormalized 
operator [0+]_r is written as 

[0 + ] R =(l + ^-m)0 + . (6.4) 

Here we only include the contributions of one-particle irreducible diagrams and ignore 
the field renormalization; the corresponding expression for H is obtained from Eq. (14.111) 
by throwing away the term 2a q , i.e. H = — 2(t b ® t b )Ti,. Below we also use the notation 

[oy R =[o\ R -o 

for the divergent (one-loop) contribution ~ 1/e so that e.g. = a s /47reHC + . 

Similar, for the "plus-minus" operators 0\ 2 we define 

\OX R { Zl ,z 2 ) =^{p^*0*](*i,*2) + W^O f ]{z u z 2 )] (6.5) 

and our task will be to find explicit expressions for the 2 -> 3 kernels H 1 ^ and H 2 ^ / . 
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6.1 transformation 

An equation for the three-particle counterterms and H 2 ^ can be obtained by the 

application of the transverse derivative — fi a d a a^ a to the renormalized operator in 
Eq. (16.41) . One obtains 

d^[0+) R = [d^+ ® iP + ] R + [tP + ® d^ + ] R . (6.6) 

We remind that a light-ray field is defined including the gauge link ip±(z) = [0, z]ip±(z), 
cf. (I3.ip . so that 

d 

d^+(z) = a^- x — [y, zn + y]ip+(zn + y)\ y=0 . (6.7) 

It is convenient to impose the light-cone gauge condition A + = on the background 
field, the same as for the quantum field a + = 0. In this case [y, nz + y] = 1 and 
^^\[y^ nz + 2/] = 0, so that the gauge links can be ignored throughout this calculation. 

Making use of the Fierz identity for Weyl spinors, (ab)(cd) = (ac)(bd) — (ad)(bc), one 
can rewrite d^x^X as follows 

d^X =[D^ + T + iff^A^i = 2d + ^t + igt\ v A b ^ - 0*A)[ASV]< , (6.8) 

where we also used that \D\ = XDX = 2d + thanks to the gauge condition A + = 0. The 
"plus" derivative can further be reduced to the derivative over the light-cone coordinate: 
d+ip-(z) = n a (d/dy a )ip-(zn + y)\ y =o = (d/dz)i/)-(zn). Thus one obtains, e.g. for the 
first term in Eq. ( 16. 6ft 

[d^+izi) ® ^+{z2)\r = <g> ?p + (z 2 )} R + ig(t b ® I^A^z^+izt) <g> ^+{z 2 )] R 

-(v\)[\DiP( Zl )®il> + (z 2 )) R . (6.9) 

Here and below we do not show color indices for the quarks and adopt a shorthand 
notation (t b ® <g> ip + ) = (4 ® 4 fc ')(V'+ ® etc. 

Combining Eqs. (16.91) . (16. 6p and (16.41) one derives for the divergent contribution ~ 1/e: 

(^1, 22)]iJ - 

= ^ A [0 + (z 1; z 2 )]' fi + (M){[A^(^i) ® lM*a)]'ji + IMzi) ® A^(z 2 )]' H } 

- ® /)[^(^)^ + (^) ® iM**)]'* - ® t b )[^ + (^i) ® ^ 2 )^(z 2 )]^6.10) 

The l.h.s. of the above equation contains the operators we are interested in. The r.h.s. 
is a bit messy so that we still have to do some rewriting. Though we are interested in 
contributions of three-particle operators only, for completeness we will keep trace of all 
(singular) terms. 

Let us start with the second term on the r.h.s. of Eq. (16.101) which contains the EOM, 
\Dip(zi). It can be rewritten as follows 

(jUA)[AZ>0(2i) ® i/j + (z2)]r = (fM\)[XDip(z 1 )] R ® ip+(z 2 ) (6.11) 
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and further 

{fiX)[XD^] R = Z_(//A)A[L% = Z_(2Z_d + ^_ - Z + D^ + ) 

= 2Z_(Z- - Z + )d + ^- + Z-Z+(n\)\DiJ) , (6.12) 

where [.D^lo means that the fields and the coupling constant which enter this expression 
are replaced by the bare ones. The factor Z + and Z_ are the field renormalization 
constants of "plus" and "minus" quark field, Eq. (I3.22p . and we also used that gA i (in 
the covariant derivative) is not renormalized. We obtain 

{ft\)[\DiP{ Zl ) ® M^)]' R = {Z-Z+ - 1^ (iJ,\)\Dip(zi) <g) ^+(^2) 

+ 2Z_(Z_ - Z + )d Zl ^{ Zl ) <g> i> + {z 2 ) . (6.13) 

The first term on the r.h.s. of ( 16 . 1 3f) is again EOM, which can safely be omitted, and 
the second term gives a contribution to the 2 — > 2 kernels but not to the 2 — > 3 kernels. 

Next, consider the contributions in the third line of Eq. (16.101) . Since the counterterms 
all have two-particle structure in one loop, one can split e.g. the first contribution, 

[AK{zx)ip + {zi) <g> i/j + (z 2 )}' r , in three terms: 

hM*l) ® V+(^)] / h4a(^) + [^+(^)4a(^)]^® + lM*l) ® [V^^l)]'*- 

The first term it given by Eq. (16.41) . The second term involves renormalization of a local 
operator. A straightforward calculation gives 

K-M' R = -7 lof^'{ " + (a^)(*W} . (6.14) 

Note that the gluon field only enters the covariant derivative, as expected in the back- 
ground field formalism. The expression in (16.141) is a one-particle counterterm which 
multiplies ij) + (z 2 ), so it contributes to the 2-^2 kernels but not to the 2 — > 3 ones. To 
handle the third term we rewrite it as 



= -z 12 ^X) f dr[M^)P + M2))' R +[M^)Al- x (z 2 )]' R , (6.15) 
Jo 

where we used the identity 

= ~ z ^) f dr p ++ (zl 2 ) . (6.16) 







The second term in (I6.15P is again a one-particle local counterterm, Eq. (16.141) . so it will 
not contribute to the 2 — > 3 kernels, and the first term can be written in terms of the 



BFLK kernel of type D in Sect. HI omitting the field renormalization. 
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To simplify the following expressions, we introduce the "averaging" operator § acting 
on a function of three variables such that 

[§ip](zi,z 2 ) = / ds (p(zi, z 2 , s) = zu / dr (p(z x , z 2 , z\ 2 ) • (6.17) 

Jz 2 JO 



In this notation 



a-.. 



-(^A)-p[SH 23 Oj](^,z 2 ) + ... 

47T6 J 



(6.18) 



where the two-particle operator H 2 3 is given in Eq. ( 14.171) omitting the terms in a q ,a g . 
The subscripts in H 2 3 specify that the operator acts on the coordinates (and color indices) 
of the second and the third field in the operator Of(zi, z 2 , Zs). The ellipses stand for the 
contributions of local two-particle counterterms, Eq. (16.141) . 

Adding the second term ~ [ip + (z 1 ) (g> ^ + (z 2 )A ^(z 2 )]' R one obtains for the full contri- 
bution in the last line in Eq. (16.101) : 



- (t fe ® J) (//A) [S i 23 Oj] , z 2 ) + (/ ig) t b ) (//A) [S 1 13 O b ](z 1 ,z 2 )} + ...(6.19) 



Last but not least, we have to deal with the first term on the r.h.s. of Eq. (I6.10p . 
The explicit expression for H + is given by an integral over the light-cone positions of 
the product of fields ip+ ® V'+j see Eqs. ( 14.111) and ( 14.11) . The transverse derivative d„x 
of this equation can be worked out with the help of Eq. (16. 9ft : 



d- x [MO + ](z 1 ,z 2 )=2 



m [d 1 o 1 + d 2 o 2 



,Zl,Z 2 ) + 



+ igU 



(t b ® I)A b -Jzi) + (I ® tVlxte) z 2 ) + EOM.(6.20) 



Finally, collecting everything and using the definitions of the kernels in Eq. (16.51) we 
obtain (omitting EOM contributions) 



Y,di[W^ k O k ]( Zl ,z 2 ) 

i,k=l 
2 

J2d k [^ f O f ]( Zl ,z 2 ) 



J2[(W + 2t b ®t b )d k O k ]( Zl ,z 2 ), 
=i 

ig\[(MA-Am)0 + }(z l7 z 2 ) 



(6.21) 



k=l 
1 



k=l 



+ (//A) [§(Jt b ® J)H 23 - (/ ® t")H 13 J C^j 22) I (6.22) 

for the contributions of two-particle and three-particle operators, respectively. Here we 
introduced the operator A which acts as 

(6.23) 



[A^(z l5 z 2 ) = (4 ® I^a(^i) + ® ^VKa(^) ¥>' j *0 



35 



This operator depends explicitly on the gluon field A. The commutator [H, A] in 
Eq. (I6.22p can, however, be rewritten in terms of the field strength tensor / + as a 
consequence of the identity [(t a <g> 7) + (I <g> t a ), (t b ® t b )] = which can easily be verified. 
Thanks to this identity the commutator [H, A] would vanish if the gluon fields in the 
first and the second term in Eq. (16.231) were taken at the same space-time point. Hence 
[H, A] ~ A b ^(z 1 ) - A\(z 2 ) and this difference can be rewritten in terms of / + using 
Eq. fl6TT6|) . ^ 

Taking into account the explicit expressions for the BFLK kernels H, H 23 , H 13 one 
can bring Eq. (16.221) to the following (final) form: 

2 3 

J2[dk^ f O f ]( Zl , z 2 ) = g(ji\) Ci[(SWi - %)O f \(zx, z 2 ) , (6.24) 

fc=l i=l 

where Cj are the color structures: 

d = f bcd (t b ® t c ) , C 2 = i(t b ®t d t b ), C 3 = -t(t d t b ®t b ), (6.25) 

the operator § is defined in Eq. ( I6.17|) . Wj are given in terms of the invariant kernels asQ 

yvi = n 23 + n 23 -n 12 -2(nt 3 + ni 3 ), w 2 = 2h 23} w 3 = 2ht 3 (6.26) 

and 

T 1 = Vi3 + V 23 , T 2 = V 23 , T 2 = Vi 3 , (6.27) 



with 



[Vi3<p](zi, z 2 ) =z 12 da d(3-(p(z" 2 ,z 2 
Jo i« " 



) ^21 J ) 



[V 2 3<f](zi, z 2 ) =zia / da dp- (pfa, zgi, z^) • (6.28) 

Jo Ja " 

It remains to solve this equation and find the explicit expressions for the 2 — > 3 kernels 
6.2 SL(2,M) decomposition 

Taking matrix elements of the operator Of over suitable states one can view Eq. (16.241) 
as an operator identity 

2 3 

£ d k U k ^ f = g(jiX) - %) (6.29) 

k=l i=l 



' Recall that the kernels depend implicitly on the conformal spins of the fields they act on. 
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on the space of functions ip^ a (zi,Z2,Zs) i— > f^izi, z<i). We will show that the kernels 
H 2 ^ are completely determined by this equation. To this end we have to examine the 
properties of Eq. ( 16.291) under the SL(2, M.) conformal transformations. 

The conformal invariance of one-loop QCD evolution equations implies that the ker- 
nels of interest, H 1 ^ and H 2 ^, are SL(2, IR)-invariant operators mapping 

H w : T 1 <g> T 1 <g> T 3/2 ^ T 1/2 ® T 1 , 

H 2 ^ : T 1 ® T 1 ® T 3/2 T 1 <g> T 1/2 . (6.30) 

The precise statement is that the operators H. k ^f intertwine SL(2, R) generators in the 
corresponding representations: Let 

£(ii,j 2 ) = ^ + ^ 2 ^ sg^) = + <?A + 5|3 ( 6 31 ) 

be the generators acting on the two-particle and three-particle states with the given 
conformal spins, respectively. For example, the first equation in f 1 6 . 3 j) means that 

H W£(l,l,3/2) = £(1/2,1) hW {6 32) 

and similar for the second operator. 

In turn, for the operators W;, S, % appearing on the r.h.s. of Eq. (16.291) one obtains 

W l :T 1 ®T 1 ® T 3/2 (-»• T 1 <g> T 1 <g> T 3/2 , 
§ : T 1 <g> T 1 <g> T 1 ^ T 1 <g> T 1 , 

T^.T 1 &T 1 &T 1 ®T l , (6.33) 

where the first equation follows from the definition of Wj in terms of invariant kernels 
introduced in Sect. 4, and the last two can be checked by explicit calculation. 

_____ ^ 2. 3/2") 

Let us apply Eq. (16.291) to the step-up operator S^' ■ On the r.h.s., taking into 
account that S^' = S^fs' + £3 one finds 

(§W, - %)Sti' im = S__ (1,1) (SWi - 7") + (S* 3 Wi - 7^ 3 ) • (6.34) 
Similarly for the l.h.s., using S 1 ^ = + Zi an d 9_T° = T 1 ^ one derives 

it=i fc=i fc=i 

Thanks to Eq. (167291) the terms in IjCT} . (16735]) multiplying S' 1 + 2 (1 ' 1) are equal; hence the 
remaining terms have to be equal as well: 

2 3 

dk z k M k ^ f = g(fi\) di^sWi - %z 3 ) . (6.36) 

k=l i=l 
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The two equations in (I6.29p and (16.361) can be rewritten equivalently as 

3 3 

8 lZl2 U 1 ^ + H 2 ^ =g(ji\) J2 Ci(S (z 3 - z 2 )W t - %z 3 + z 2 T t ) = J2 CiA , 

1=1 i=l 
3 3 

d 2 z 21 U 2 - f + H 1 ^ / =g(fi\) c *($ (*3 - *i)W, - %z 3 + zi%) = °i B i ■ ( 6 - 37 ) 



i=l i=l 



These equations are already SL(2,M) invariant and correspond to the mappings T 1 <g> 
T 1 (g>T 3 / 2 h+ T 1 OT 1 / 2 and T 1 <g> T 1 <g> T 3 / 2 ^ T 1 / 2 ® T 1 , respectively. The expressions on 
the r.h.s. can be simplified. After a straightforward calculation one derives the following 
expressions 



[Aicplfa, Z2) = z\ 2 ( / df3f3(p(z 1 ,z 2 ,zf 2 ) - da dp p f(z 1 , z% v z{ 2 ) 
\Jo Jo Jo 

[A2(p]{zi t z 2 ) = z\ 2 da dp — <p(z lt z% v zf 2 ) , 

Jo Ja « 



1 



a 



[A 3 <p](zi, z 2 ) = z l2 I da I d p—{a ~ P) V( z i2, z 2, 4i) 



1 



\B\<p]{z\, Z2) = ~z l2 / d(3(3 (p(z\, z 2 , z' 21 ) - da df3 (3 <p{z" 2 , z 2 , %) 
f 1 f 1 a 

[B 2 (p]{z!, z 2 ) = -z\ 2 da d{3— (a - p) (p(z 21 , z 2 , z{ 2 ) , 
Jo Ja '> 

[B 3 (p](zi,z 2 ) = -z\ 2 / da / dp — tpfa, z™ 2 , z%i) ■ ( 6 - 38 ) 
Jo Ja a 

As the last step, we can prove that the two equations in (16.371) or, equivalently, 
Eqs. (I6.29p and (I6.36P determine the kernels H fc ^ uniquely. Using the first equation in 
(16.371) to eliminate H 2 ^ one obtains for 

3 

(d 2 z 21 d lZl2 - l)U^ f = g(fiX) °i (d2Z 21 Ai - Bi) . (6.39) 



i=i 



The differential operator on the l.h.s. of this equation is nothing else as the Casimir 
operator for the conformal spins (ji, j 2 ) = (1/2, 1): [d 2 z 2 \d\Z\ 2 — 1) = S\ 2 — 3/4. This 
operator has zero modes (Sf 2 — 3/4) ^ = 0. However, functions of the form ip = M l ~*fip 
do not belong to this subspace Therefore, this operator can be inverted and 



* Indeed, Eq. (|6 .32[) implies that H 1 ^ maps the eigenspace of the Casimir operator Sf 23 onto the 
eigenspace (with the same eigenvalue) of the operator Sf 2 . Since, the minimal eigenvalue of (<S' 2 ) 12 ' 3 1 ' 3 ^ 2 ' ) 
is 35/4, one concludes that if ip = U^l'tp then ((S' 2 )^ /2 ' 1) - 3/4)</> ^ 0. 
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restored as 

3 

H 1 ^ = g(ji\) CM (d 2 Z2iA - Bi) , (6.40) 
i=i 

where the IR is the inverse operator to (Sf 2 — 3/4) on the subspace that excludes zero 
modes, {Sf 2 - 3/4)^ ^ 0: 



{Ri}}( Zl ,z 2 ) = 




(6.41) 



In practice, this last transformation is not necessary as the solution to Eqs. (16.371) can 
easily be guessed. The final result is presented below in Eq. (17. 7p in Sect. [3 It coincides 
with the corresponding expression in Ref. obtained by a direct calculation. 

As the final remark, we want to note that the case of quark-quark operators with 
the same chirality considered here proves to be the most complicated in this respect. In 
most other cases solving the similar pair of equations is straightforward since, as it turns 
out, one (or both) of the 2 — > 3 kernels enters without a derivative. 



6.3 M MM transformation 

Equation (I6.36P that we derived using conformal SX(2,R) transformations can also be 
obtained in a different way, by applying the generator of Lorentz rotations M MAt to the 
quasipartonic operator + . This alternative derivation may be interesting as it shows 
that Lorentz symmetry alone is sufficient to restore the "nondiagonal" 2^3 operator 
mixing, so in what follows we outline the main steps. 

The general strategy is similar to the case of the transformation considered in 
Sect. 16.11 Using the definition in flA.ll) one obtains 

1 d 

i[M^,ip + (z)) = -(fi\)5 M ip+( z ) , <W = 2 zd n\ + V-qT ■ ( 6 - 42 ) 

The first difference to the P^^ case is that contributions from the rotation of the gauge 
links in the definitions of light-ray operators are nonzero and have to be taken into 
account. Explicit calculation yields 

5 M [0, z]il>+(z) = (zd z + l)if)-(z) - % ^rz 2 {fi\) [ duuf ++ {uz)ij + { Z ) . (6.43) 

^ Jo 

A bigger problem is that 5m is not a symmetry transformation of the QCD Lagrangian 
because of the gauge fixing term 

i j rf 4 x £„ = ~hj d * x ^^(x)f . (6.44) 

Applying the Lorentz rotation to the path integral, one obtains for a generic com- 
posite operator O 

[SmO}' r = 6 m [0]' r +[(6 m S)0]' r (6.45) 
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i.e. there is an extra term due to the variation of the QCD action § 



5mS = — - J d 4 x a + (x)5Md+(x) . (6.46) 
In our case one obtains 

5 m O + ( Zi , z 2 ) = (z 1 d 1 + z 2 ) + {z 2 d 2 + l)0 2 { Zl , z 2 )- 

Tdr{zl(t b ® l)0){z u z 2 , tz x ) + z\(l ® t b )O b f (z u z 2 , tz 2 )} + EOM 

so that 

(sift + l)[0 1 ]' R (z 1 ,z 2 ) + (z 2 d 2 + l)[0 2 ]' R (z 1 ,z 2 ) = 

rdr[zl{t b <g> I)[O b f ]' R { Zl , z 2 , tz x ) + z\[I ® * 6 )[0j]' B (*i, *2, tz 2 )} 
+ 5 M [<D + ]' R (z u z 2 ) - [0 + (5 M S)]' R + ... (6.47) 



19 , 
2 



cf. Eq. (16.101) . Note the l.h.s. of this equation is the same as in Eq. (16.361) . The r.h.s. 
is similar to the r.h.s of Eq. (16.101) except for the additional term [0 + (5mS)]' r . Let us 
consider this contribution more closely. 

The gluon propagator corresponding to (16.441) is in general (in vector notation) 



G^ v {k) - — 



-9^u + 



kfjiy + n u k 



kn 



(6.48) 



and the light-cone gauge corresponds to taking limits n 2 — > and £ — > 0. 

The Lorentz rotation of the action (I6.46P is formally 0(1/ '£). However, since n^G^ = 
—i^k u /{kn) (16.481) . each insertion of a+(x) is effectively 0(£), so that terms in 8m S 
containing a+ more than once can safely be neglected. Hence it is sufficient to keep the 
first term only in the expression for the transformed field 



S M a+{x) 



1 / \ 1 
:Sa( x ) + 



2"^-' ' 2(ji\) 
and consider zero-momentum insertions of the operator 

1 



(nxdp,)a + (x) 



>(ymrtr = — a + (x)a^(x). 



(6.49) 



(6.50) 



A detailed analysis shows that the effect of such insertions reduces to the following 
transformation of the light-ray fields: 



(5 M S)[0,z)i/> + (z) 



du a„x(uz) 



a 



i/> + (z). (6.51) 



§ Strictly speaking there are also terms due to breaking of Lorentz symmetry of the counterterms 
to Cqcd itself, i.e. field and coupling renormalization. Such corrections generally correspond to one- 
particle-reducible self-energy insertions and are of no relevance here. 
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Figure 2: 6m S insertions, see text 



This equation is illustrated in Fig. [2j The first term on the r.h.s. comes from the gauge 

link, Fig. [2k. Note that a + (x)a + (y) = —i^5^{x — y) so that the insertion of 5mS simply 
replaces a + {uz) by the physical transverse gluon field a^\{uz) at the same position. The 
second term on the r.h.s. of ( 16.511) corresponds to the first contribution in Fig [2b, which 
originates from the contraction of the quark propagator. There is also another term that 
involves the classical background field A \ and arises when one completes the covariant 
derivative to obtain EOM. Such contributions cancel, however, in the sum of all one-loop 
Feynman diagrams so we do not show this term in Eq. (16.511) explicitly. We have checked 
that all other insertions of 5mS (and also the same insertions as in Fig. [2] but with the 
contraction of the transverse field a^x) m the effective vertex instead of a + (x)) cancel 
as well. The structure of these cancellations strongly suggests that they can be put in a 
form of a certain Ward identity which we did not work out in the operator form, however. 

Using (I6.16P the contributions ~ in Eq. (16.511) can be rewritten in terms of the 
field strength tensor / ++ so that the last term in Eq. (16.471) . [0 + (SmS)]' r , reduces to the 
sum of contributions of the corresponding BFLK kernels. The rest of the calculation is 
straightforward and follows closely the calculation described in Sect. 16.11 so we skip the 
details. The result reproduces Eq. (I6.36[) . 



7 Results for the 2^3 Kernels 

In this Section we present a complete list of the 2 — > 3 kernels. Altogether, there exist 
72 different twist-3 pairs of the primary fields, {ip- <g> ip + , ip- ® x+i 4'+ ® • • •}• 

Many of them are related to each other by parity P and charge conjugation C symmetries, 
however, so that there are only 16 independent pairs. Eight of them can be chosen as the 
descendants of the quasipartonic X + ® X + operators with both fields having the same 
chirality and the other eight as the descendants of quasipartonic operators containing 
primary fields of opposite chirality. In each case, the renormalized operator [X]r(zi, z 2 ) 
in the light-cone gauge can be written to the one-loop accuracy as 

[X] R (z l7 z 2 ) = X B ( Zl ,z 2 ) + ^L[M^X](zi,z 2 ) + ^[M^Y]( Zl ,z 2 ) . (7.1) 

47T6 47T6 

The results for the 2^2 Hamiltonians, M^ 2 \ are collected in Sect. [5j The ex- 
pressions for the 16 independent three-particle counterterms [E[( 2 ^ 3 )y](z 1 , z 2 ), where 
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Y(zi, z 2 , z 3 ) = X+\z\) ® X+\z 2 ) <8> X+\zs) is a quasipartonic operator with proper 
quantum numbers, will be given in what follows. 

The representation in (17.11) is somewhat schematic since several three-particle op- 
erators Y±, Y 2 ,...can contribute to the r.h.s. and all have to be taken into account, 
which makes the corresponding expressions quite cumbersome. For example, the oper- 
ator / + _ g) if) + can mix with the operator / ++ (g> ® / ++ and also with a three-quark 
operator ip + ® ip + ® etc, so that, effectively, the gluon pair / ++ ® / ++ is replaced by 
the quark-antiquark pair ip + ® Such quark- ant iquark pairs always appear in special 
combinations 

J a (z u z 2 ) = i> + {z l )T a ij + (z 2 ) + x + {z 1 )T a x + (z 2 ) = q( Zl )T a 1+ q(z 2 ) , J a (z) = J a (z,z) , 
J a \z 1 ,z 2 )=Mz 1 )T a T b ^2)-X + (^)T b T a x + (z 1 ) } 

J abc ( Zl , z 2 ) = M Zl )T a T b T c M^) + x + (z2)T c T b T a X+ (zi) , (7.2) 

where the sum over flavors is implied, ip + T a ip + = s ^2 lA ip^_T a ip^_^ etc. The operators 
J a , J ab , J abc are of course related to each other. However, imposing such relations com- 
plicates the expressions considerably, so we do not attempt this. 

As mentioned above, sixteen independent 2 — > 3 RG kernels will be presented explic- 
itly. The remaining 56 = 72 — 16 ones can be restored as follows. 

First, the left-handed and right-handed quark spinors appearing in (17. ip can be inter- 
changed freely: ip± x± an d $± X± ( on both sides of Eq. (17. ip simultaneously, but 
separately for each flavor). The form of the kernels is not affected by such transformation 
so that, for example, the four kernels 

xp- ® ip + -> ip. ® ® /++ , X- ® ^+ -> X- ® ^+ ® /++ , 

^_ ® x+ -»• ^- ® X+ ® /++ , X- ® x+ -> X- ® x+ ® /++ 

are equal to each other 0. Making use of such simple substitution rules one ends up with 
36 kernels: 16 = 4 x 4 quark-quark, 16 = 2 x 8 quark-gluon and 4 gluon-gluon ones. 

The remaining 36 = 72 — 36 kernels can be obtained by applying hermitian con- 
jugation to the corresponding quantum operators. To this end, the color factors like 
i(t d t c ) ® t c should be treated as c— numbers, i.e. 

(i(t d t% t , ® t] r ^( Zl Wl(z 2 )) ] = [i{tH% v ®t^\^{z 2 ))\^{z^ 

Here in the first line all SU(N) generators are taken in the quark representation, t a = T a , 
whereas in the second line t a = (—T a ) T are the generators in the antiquark representa- 
tion, cf. Eq. (ET5|) . 

^ The annihilation-type contributions ~ Sab in Eq. (|7.43p are only present for ip^tp^, x+X- opera- 
tors and have to be discarded for x+V'-i V'+X-- 
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We remind that under hermitian conjugation 

4 = ^ ± , xl = X±, fl + = f++, /!- = /+-, = (7.3) 

and in addition one has to replace (//A) — > (A//)- One also finds easily that 

( z 2 ))t = r(z 2 , Zl ) , ( J afe (^, z 2 ))t = j ba ( Z2 , Zl ) , ( j abc ( Zl , Z2 )y = r ba (z 2 , Zl ) . 

(7.4) 

7.1 ^-^+ or ^D- + ip + 'ijj + 



X^{z 1 ,z 2 )=^_{z l )^ + {z 2 ), 

Y ijd ( Zl , z 2 , z 3 ) = g{li\)1>%{ziW + {z*)fi+{zz) 

X i i(z 1 ,z 2 ) = ±D_ + 4> i + (z 1 )ft + (z 2 ) 



The three-particle counterterm in both cases takes the form 

3 



k=l 

where are the color structures 

(Ci)% d = f h %t) 3 , , (C 2 )%, d = i(fV)«4, , (C 3 )^ = <4(^)«' • (7-6) 
The 5X(2, M) invariant kernels for the first case, X tj = V'LV'+j are given by: 

pi pa 

\Hiip\{z u z 2 ) = z\ 2 I da j dpptp(z* 2 ,z 2 ,z 2 3 1 ) , 
Jo Jo 

[HM( Zl , z 2 ) = z\ 2 f da f dfi^ V (z« 2 , z 2 , 4) . (7.7) 

Jo Jd « 

(the third color structure, C3, does not contribute). 
For the second case, X tj = | L>_ + ■?/>!,_ t/^, one obtains 

[Hiip](zi, z 2 ) = z 12 \ I dp (3 ip(z!, z 2 , z{ 2 ) + da dp (3 ip(zi, z" 2 , z 2l ) 

'0 Jo Ja 

1 r i m 

P ..t..a „ J 







- / da j dp— cp(z^, z 2 , z p 12 ) } , 



[H 2 (p](z 1: z 2 ) = z 12 f da f dp— (2 - ^) <p(zn, to, 4i) » 
Jo Ja « V «P/ 

[H 3 (p}{zi,z 2 ) = -z 12 da dp P (p(zi, z^, z{ 2 ) . (7.8) 

Jo Ja 
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7.2 \D_ + f ++ U 



X a \ Zl , z 2 ) = \D-+fl + {z x )fX + {z x ) , Y M (z x , z 2 , z 3 ) = gi^flJzJftJzjEjzs) 



The three-particle counterterm takes the form 

3 

i.TTT 



k=l 



where the color structures are the same as in Eq. (17.61) (but with all generators in 
the adjoint representation, cf. (I3.3P ). The SL(2,M) invariant kernels are: 



[HM(zi,z 2 ) = z 12 { I d/3 (3 (p(z 1} z 2 , z'( 2 ) + / da j d/3— I 2 - — ) (p(z^ 1} z 2 , zf 2 



JO J a 

1 r± 



+ J^da J_ dPapU-^tp(z u % 2 ,4 x ) 



[H 2 (p](zi, z 2 ) = z 12 I da [ d(3^-(p(z™ 2 , z 2 , z^) , 
'o J a « 

\Hzip\{z x ,z 2 ) = -z X2 \ da dp <p(z u z 2l , z{ 2 ) . (7.10) 

J0 Ja « 



7.3 ±D_ + $+/ ++ or $ + \D_ + f ++ 



X i \z 1 ,z 2 ) = \D_ + mz 1 )jlAz 2 ) 



or z 2 , z 3 ) = gfaX^iizJf'+Wft+iza) 



X ia (z l ,z 2 )=^Uz 1 )\D„ + flJz 2 ) 



The three-particle counterterm takes the form 

5 



47T€ 

fc=l 



The first three color structures Ci, C2, C3 are the same as in Eq. (17.61) (with the generators 
in the appropriate representation), and there are two new structures: 

= i{t d e't a ) lV , (c 5 )%, d = i(e't d t a ) u , , (7.12) 
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The invariant kernels for the first case, X ta = ^D- + tp l + f+ + , are given by: 



[Hi(p](zi, Z2) = ^12 1 y df3f3ip(z 1 ,z 2 ,zf 2 ) + y rf«y d/3 jj- ¥(2%!, Z2, z% 2 ) 

[W2¥>](^1, ^2) = 212 j Q da j_ d fi~~ ( 2 ~~ ^(*12> ^2, 22i) , 

[H 3 (p}(z 1 ,z 2 ) = -z 12 / da / d/3 <p(zi, zf 2 ) , 

Jo Ja a 

[H 4 (p](zi,z 2 ) = -z 12 da d[3a[3ip(z 2 ,z* 2 ,z 21 ) , 

JO Ja 
pi pa 

[H 5 ip](zi, z 2 ) = -z 12 / da d(3 a(3 ip(z 2 , z" 2 , z^) . (7.13) 
Jo Jo 

The invariant kernels for the second case, X ia = ^(zi)|Z}_ + /" + (z 2 ), are: 

[HM(z u z 2 ) = z 12 \ [ dpp<p(zi,z 2 ,z$ 2 )+ f da f dp (5 (p(z? 2 , z 2 , z^) 
I Jo Jo Jo 

[H2(p](z u z 2 ) = z 12 / da d(3 (3ip(z* 2 ,z 2 ,z 2 3 1 ) , 

JO Ja 

[H 3 <p](zi, z 2 ) = -Z12 / da dp ifiz!, z 21 , z^ 2 ) , 

Jo Ja a 

[H 4 ^}(zi,z 2 ) = -z 12 da d(3a(3ip(z 2 ,z^ 2 ,z 21 ) , 

^0 ia 

[n 5 <p](z 1 ,z 2 ) = -z 12 da dPaP<p(z2,z% 2 ,4i)- ( 7 - 14 ) 
Jo Jo 

7.4 or i J D_ + / ++ ^ + 



X° i (z 1 ,z 2 ) = f$. + (z 1 )rl>i(z 2 ), 

or Y^,^) = gWfl + {z 1 )tl)%{z 2 )fi + {z z ) 

X a \ Zl ,z 2 ) = \D_ + fl + { Zl )^ + {z 2 ) 
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The three-particle counterterm takes the form 

= ^ £(^SwW^* d . (7.15) 
fc=i 

The color structures C X ,C 2 , C 3 are specified in Eq. (17.61) (with the generators in the ap- 
propriate representation), C 4 , C 5 are given in Eq. (I7.12p . and there is one new structure: 

(C 6 )^ d = i(t a 't a t% . (7.16) 

The invariant kernels for the first case, X a% = fl,ip l _ are given by: 



[Hiip](zi, z 2 ) = -z 12 { I da J dp (3 tp(z u zf 2 , z% x ) 



+ f da I df3 [ drya-f (l- ^) tp(z? 2 , 

Jo Jo Jb V °nj 



P 7 
Zoi i Z' 



21 1 *21 1 ( i 



[H 2 tp}(z 1 ,z 2 ) = -z( 2 I da I dpi d'y —<p(z? 2 , z^, z^) 



JO 



a 



[H 3 (p]{zi, to) = z\ 2 1 jf da J dp^-ip{zt, z 21 , z{ 2 ) (7.17) 
+ C da f dp f dryer/ (2-^) ip(z? 2 ,zi,zl) } 



Jo Jo V "7 



[Ha^Uzi, z 2 ) = z\ 2 / da I dp I ^707 ( 2 - ^-M 
Jo Ja Jo V °nJ 



^(^12) z 21i z 



21 J ' 



[H 5 (f](z u z 2 ) = z\ 2 I da I dpi d'y — <p(z X2 , z^, zl x ) , 



J a 



a 



n\ nl n\ — 2 — 

[Help] (z u z 2 ) = -z\ 2 / da dp dj — (p(z X2 , z^, z^) 
Jo Ja Ja " 



(7.18) 
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For the second case, X ai = \D-+j\^\ one obtains 



[HM(*i, **) = *di f dPP<p(z u z 2 , 4) + J 1 da J 1 d(5^-{2- ^ <p(&, z 2 , z{ 2 ) 



+ [ da fdp ^ d 1 a 1 { 4 - ^) <p(z? 2 , 4, 4i) 

JO JO J 13 



J da J d/3 —ip(z" 2 , z 2 , 4) 
+ £ da £ dp £ rf 7 «7 (2 + ^ </?(^ 2 , 4, 
[H 3 (p}(zi, z 2 ) = -z 12 J da J dp J d'yory^A-^j <p(zf 2 , 4, ^21) , 
[ft^K^i, 22) = -*i2 jf da ^ d/5 ^7 «7 ^ - ^0 ip(z^ 2 , 4, ^21) > 



c?7 «7 ( 2 + — j <p(z? 2 , 4, 4i) , 
cry 



\H b (p\{z 1 ,z 2 ) = -z l2 I da I dp I 

Jo J a J a 

[Hm\(z x ,z 2 )=z 12 J\a J\p f^ar/ (2 + ^ ^,4,^). (7-19) 



7.5 


/+_V>+ or 












fl_( Zl W + (z 2 ) 


Y aid ( Zl , z 2 , z 3 ) 


= g^\)fl + ( Zl )^ + (z 2 )fi + (z 3 ) 






or 


J m (z u z 2 ) 


= g(fi\)J a (z 1 )^ + (z 2 ) 




X m (z u z 2 ) = 


fl + { Zl )i,l{z 2 ) 


Z abi (z 1: z 2 ,z 3 ) 


= g^\)J a \z ll z 2 )^ + (z 3 ) 



The three-particle counterterm for the first case, X ai = fl_vp 1 ,, takes the form 



Itvc ' 



[X a % = ( E( C7 *)^ W * yB ' < ' d + ^1^' + i^Z aU ' ) , (7.20) 



vfc=l 
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where the color structures G\,C<i are specified in Eq. (17. 6ft . The invariant kernels are 

pl pa 

[HM(zi, z 2 ) = z\ 2 \ da dp a(3 <p(z? 2 , z 2 , z^) , 
Jo Jo 

[H2<p]{zi, zi) = z\ 2 \ da dp (p(zi 2 , z 2 , z^) , 

Jo J a " 

[HitpKz!, z 2 ) = z 12 daa 2 (p(z" 2 ,z 2 ) , 
Jo 

[H 2V }( Zl ,z 2 ) = -z 12 J daj^ dpi (7.21) 
For the second operator, X ai = f+ + tp l _, we obtain 

[x a % = £^ (j2(c k ):%, d n k Y^' d + (to'ehm^' + ti,n 2 z b ^ , (7.22) 

where the color structures C k are defined in Eqs. fl7.6l) . fl7.12p . The invariant kernels are 

pl pa 

[Hi<p]{zi, z 2 ) = -z\ 2 I da j dp P <p(zi, z 21 , z{ 2 ) , 
Jo Jo 

[H 3 <p](zi, z 2 ) = z\ 2 I da I dp — (p(z u z 2l , z{ 2 ) , 



J a 

1 



a 



pi pa 

[Ha<p]{zi, z 2 ) = -z\ 2 \ da dp a/3 ip(z" 2 , z h z^) , 
Jo Jo 

pl pa 

[H 5 <f}(zi, z 2 ) = -z\ 2 \ da dp aP ip(z^, z u z{ 2 ) , (7.23) 
Jo Jo 



(the second color structure does not contribute) and 

[HiLp]{z 1 ,z 2 ) = z 12 I da aa ip{z^ 2 , zx) 
Jo 



pl pa 

[H2f](zi, z 2 ) = Z12 / da dp <p(z? 2 , z 1: z^) . (7.24) 
Jo Jo 



7.6 ^ + / + _ or \D_ + $ + f. 



X ia {z 1 ,z 2 )=i>\{z 1 )fl_{z 2 ) Y™ d (z u z 2 ,z 3 ) =g( fJ ,\)ft + (z 1 )ft + (z 2 )f* + (z 3 ) 
or J ia ( Zl , z 2 ) = g{fxX)fi{ Zl ) J a (z 2 ) 

X ia ( Zl ,z 2 ) = \D_ + i,X{z x )f^{z 2 ) Z iab (z u z 2 ,z 3 ) = g^\)r + { Zl ) ® J ab (z 2 ,z 3 ) 
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The three-particle counterterm for the first case, X ia = ip+fl-, takes the form 



fe=l 



(7.25) 



with the color structures Ck as defined in Eqs. (17.61) . ( 17TT21) . (17. 161) . The invariant kernels 
are 

[Hi(p](zi,z 2 ) = -zfJj da J dp ap <p(zi, z% ± , zf 2 ) 
+ [da f d(3 f a di{p<y + Pi) 

JO JO J/3 

/*1 pol rl 

[H2<p](z 1 ,z 2 ) = -2$ 2 da dp dry (^7 + Pi) <p{z?2, 4i,4i) 

JO JO Ja 

dp — (p{z!, z 21 ,z^ 2 ) 



[H-M(zi,z 2 ) = z\ 2 \ J da J 



a 

+ 2 f da ^ dp I d 1 p 1 <p{zi2,zL*li) ) • 
Jo Jo Jo 

[H4<p](z 1 ,z 2 ) = -2$ 2 da dp (jfy + Pi) (p^, , 

Jo Ja J p 

-i pi rP 

djPj <p(z" 2 , z$ u z1 x ) , 



[H s <p](zi, Z2) = -2^ 2 f da f dp f 

Jo Ja Jo 

/*1 rl i*a 

[H 6 <p](z 1 ,z 2 )=2z 2 2 da dp d 1 p 1 i P (z^z^zl) (7.26) 

Jo J a Jo 

and 

[HM(z 1 ,z 2 ) = -z 12 \ [ dpp 2 V ( Zl ,z^) + 2 I da I dppp(p(z^,4 



21/ ( 1 



JO JO 

1 r l 



[H 2 (p](zi, z 2 ) = -2zi 2 / da dppP(p(zu,Z2i)i 

Jo J a 

[H 3 ip)(zi, z 2 ) = 2i 2 |^ da J dp jj<p(zi, z° 2 , zjj + da J dp J djip(z? 2 , z^, zj^j . 

(7.27) 
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For the second operator, X ia = ^D_ + ?p l + f+_, we obtain 



,k=i 



+ t\ v U z Z l ab + t\ v H±Z l ba j , (7.28) 
with the same color structures C k as in the first case. The invariant kernels are 

1 r l r l 32 



and 



[Hitp](zi, Z2) = zi 2 \ J dp p (p(zi, z 2 , ^12) + y da y d/3— ^(^, z 2 , zf 2 ) 
+ 2 /'da f d/3 f a d>y(2p>y + P*() ( p(z? 2 ,4 1 ,zZ 1 ) 

JO JO J 8 

[H 2 <p\(z u z 2 ) = 212 jy£ rf/3 ^( 2_ ^f) ^"2,^2, 4) 

+ 2 /'da f d/3 / 1 d 7 (2^ 7 + /37)^(^2,4^2i) 

JO JO Ja 



[Hs(p](zi,Z2) = -6zi2 f da f dpi d'y P7 (p(z" 2 , 4, z^) , 

Jo Jo Jo 

[HM(z 1 ,z 2 ) = 2z 12 [ da! dpi d 1 {2p 1 + p^) v (z<? 2 , z^, zl) , 

Jo Ja J 8 

nl nl n8 

[H 5 <p] (zi ,z 2 ) = 6^12 / da / d/3 / d<y(3~f ip(z" 2 , 4, z^) , 

Jo Ja Jo 



rl nl /*a 

[He<p](z 1 ,z 2 ) = -6z 12 da d(3 d-yPytp^, , (7.29) 

JO Ja JO 

[R 1 <p](z 1 ,z 2 ) = 6 f da f dPW<p{%2,4i), 
Jo Jo 

[HM(z 1 ,z 2 ) = 6 [ da f dp 00 

JO Ja 

pi pa pa 

[H 3 <p](z 1 ,z 2 ) = -2 / da / d/3 / d 7 </?(-^2,4,4), 

JO JO J/3 

[H4f](zi, z 2 ) = - / da </?(zi, 2^,2:2) • (7.30) 
Jo 



50 



7.7 /+_/++ 



Y abd (zi, Z 2 , Z 3 ) = ^A)/^*!)/^) j?+(*3) 

yab/', „ n _ fa ( v \tb r v \ G ab (z 1 ,z 2 ) = g(n\)J a (z 1 )fl + (z 2 ) 

{ Zl ,Z 2) -U_[ Zl )} ++ [Z 2 ) E a b(zuZ2}Z3) =g(j i \)fi + ( Zl )J°<*( Z2 , Z3 ) 

Z abc ( Zll z 2 , z 3 ) = g{iiX)J a \z u z 2 )flAz 3 ) 



The three-particle counterterm takes the form 



\x ab ]' R = ^- ( J2( c ^ d n k Y a,b ' d + (c, ® ti v )HxG a ' h ' + n 2 E ab 
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fe=l 



+ t c w HsZ Mb '+t c ab/ n,Z M 'j, (7.31) 
with the color structures C±, C 2 as defined in Eq. (17. 6p . The invariant kernels are in this 



case 



and 



rl pa 

[H 1 ip](z 1 ,z 2 ) = z\ 2 \ da dp acfl ip(z" 2 , z 2 , z% x ) , 
Jo Jo 

[H 2 ip]{ Zl , z 2 ) = z\ 2 \ da dp (p(z&, z 2 , z$\) > ( 7 - 32 ) 

Jo Ja a 

[Hi^p](zi,z 2 ) = z 12 daa 2 ip(z° 2 , z 2 ) , 
Jo 

\H 2 y\{z x , z 2 ) = z 12 da d(3 a (p(z X2} z 21} z 2 ) , 

Jo J a 

[H 3 (p]{z 1 ,z 2 ) = -z 12 [ da f dp ^(p{z$ u z{ 2 , z 2 ) , 

Jo J a P 

pi pa 

[H 4 <p}(zi,z 2 ) = -z 12 da dp (3 ip(z 21 , z 2 , z{ 2 ) . (7.33) 
Jo Jo 



7.8 or \D. + f ++ f ++ 



Y M ( Zl , z 2 , z 3 ) = g^\)fl + ( Zl )f b ++ (z 2 )f* + (z 3 ) 
X ab ( Zll z 2 ) = fl + { Zl )f b + _{z 2 ) G ab ( Zl , z 2 ) = sGuA) J\z 2 ) 

or 

X ab (z u z 2 ) = \D_ + fl + ( Zl )f b ++ (z 2 ) E ab ( Zl , z 2 , z 3 ) = g{ji\)fl+(zi) J bca (z 2 , z 3 ) 

Z abc (z u z 2 , z 3 ) = g^X)fl + (z 3 )J bc (z 2 , z 3 ) 
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The three-particle counterterm for the first case, X ab = fl + fi__, takes the form 



fc=i 



+ ft 3 £ a& + C^4^ a bc + t c ba ,H 5 Z a ca j . (7.34) 
The color structures Ci,C 2 , C 3 are defined in Eq. (17. 6p and C 4 , C 5 , C 6 are given by 

(C-Oa'&'d = / de ^L *a&' > (CjOa'b'd = ^)ba> ^ a b< > (^a^'d = {^)ab> ■ (7.35) 

Note that these structures are different from those in Eqs. (I7.12p and (I7.16p . The invari- 
ant kernels are 

*1 pa 

\li v ^\(-. y . z 2 ) = -z\ 2 { j da J d(3 a/3 <p(zi, z 21 , z( 2 ) 

-1 rot ca 



+ l**J° w f d ^ *fc ( 2 - + 4, 4) 



'/3 

[W 2 ¥j](«i, 22) = -z( 2 I da f dp f d-y af3j ( 2 + ^ ) <^( z « ? ^ ^) 



o v/o -'a 

1 /■! ~,2 



a/3 



+ jf * da £ dp J' d 1 7 ^2 ap (2 - ^ + a/3 ( 2 ~ ) 4. 4) 

[n 4 (f](z u z 2 ) = -3zf 2 [ da f dpi dj <p(z? 2 , z^, Z&) , 

JO « 

[W5¥>] 22) = 34 jT da jf 1 d/3 ^° d 7 (2 - |0 p(z? 2 , 4, ^1) , 

[Hfj<p](z 1 ,z 2 ) = -3zi 2 / da / d/9 / cfy ^(212,4^21) ( 7 - 36 ) 

Jo J a Jb " 
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and 

• 1 pa 



21) f ' 



[H 1 if}(z 1 ,z 2 ) = -z 12 y^ daa 2 <p(z u z« 2 ) + 2^ da £ d/3a/3p (2 + <p(z« 2 , z{, 

pi pot 

[H 2 <p](zi, z 2 ) = -6z 12 da d(3 at/3/3 (p(z 21 , z{ 2 ) 

Jo Jo 

/*1 rck /*1 

[K^p){zi, z 2 ) — 2zu / da d(3 I d^aip( 
Jo Jo Ja 

[H^Kz-l, z 2 ) = da J_d(3^<p(zi,Zi 2 ,z 2 i) + 2 J da J^d(3 J d~f a(p(z° 2 , zj, z^) 

[H 5 ip](zi, z 2 ) = -2z 12 I da j d/3 \ d^a <p{z* 2 , z^, z^) . (7.37) 

Jo Ja Ja 



For the second operator, X ab = |-D-+/" + /+ + , we obtain 

= £^ ( f2( c ^ d n k Y a ' b ' d + (c< ® tt v ) (h.g^ + H 2 G b ' a '] 

716 \ k=l 

+ n 3 E ab + n A E ba + t c aa , (n 5 z a ' cb + n 6 z a ' bc ^ + t c ba ,n 7 z a ' ca ^j . (7.38) 

The color structures in this expression, Ci,C 2 , C 3 and C 4 , C 5 , C 6 , are defined in Eqs. (17.61) 
and (I7.35p . respectively. The invariant kernels are given by 

[HM(z u z 2 ) = z 12 < f d/3/3(p(z 1 ,z 2 ,z( 2 )+ [ da f d(3 f 2 - ^ ) (p(z^, z 2 , z{ 2 ) 

Jo Jo Jot 



f3 \ a/3 



2 da f d/3 l' X d 1 ah (4 + ^-^1 + 2^) ^z^z^zl)}. 



JO J /3 



a/3 a^j /?7 



[Ti 2 <p](zi, z 2 ) = z 12 { j da J d/3— (p(z" 2 , z 2 , z^) 



+2 L da l? p L di a ( 2 + 1) + ^ ( 2 + I) ) zL zl) ] ' 

[Ha<p](z 1 ,z 2 ) = -6z w f da f d(3 I d-yafa (2 + ^ - ^] V {z^z^z\) , (7.39) 

Jo Jo Jo V a/3 aryj 
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[7U(p] (zi, z 2 ) = 6z 12 da J dp J fa afa ^1 + (p(z^ 2 , z^j) , 
[H 5 (p}(zi, z 2 ) = -6z 12 J da j dp J fa afa (s - <p(z? 2 , z^, z^) , 
[7U<p\(z 1 ,z 2 )= 6z 12 j^da j^dp j dyafa^l + ^j^z^z^zl). (7.40) 



and 



[HM( Zl , z 2 ) = 6 jf 1 da £ dp aPP (3 + |0 <^ 2 , 4) , 

[H 2 <p\(z u z2) = 24 I da f dPaPP V (z^,zf 2 ), 
Jo Jo 



pi. rex M 

[H 3 (p](zi, Z2) = -6 / da dp fa a ip(z? 2 , z$ x , z^) , 

JO JO Ja 

[H±ip\{z u z 2 ) = - I daa(p(z? 2 ,Z!,z 2 ) , 
Jo 

[H 5 (p](zi, Z2) = - da(p(z 1: z° 2 ,z 2 ) , (7.41) 
Jo 

/*1 /'CI 

[We^K^i,^) = -6 / da / tZ/3 / ^70^(2^, 4^21) , 
Jo Jo J/3 



'0 JO J 13 

-1 r l r l 



[H7<f](zi,z 2 ) = 6 / da / d/3 / faa<p(zf 2 , z^, z^) . (7.42) 

JO ia J/3 

7.9 ^ + ^_ or iL>_ + '0 + '0 + 



-^AB^l, 22) = V'+^l)^- 5 ^) ^4s(^l, ^2, Z 3 ) = g{^W+ A {z 1 )i) 3 f(z 2 )fl + (z z ) 

or z 2 ) = g(fi\)J a (z 1 ) f b ++ {z 2 ) 

Xj B (z u z 2 ) = \D. + ^ A { Zl Wf{z 2 ) G abc (z 1 , z 2 ) = 2(M)/^0/* + (^)/; + (*3) 



Note that in this case we display the flavor indices A, B of the quark fields. 

The three-particle counterterms for the both operators, X l j B = ^ A ip 3 2 B and X % i B = 
^D_ + iJj+ a 4>+ B ', have the same structure: 



[ X ab\'r 



y^X C k)\ 3 /jld UkY l AB d 
k=i 



+ S AB ( i(t a t b ) ji H 1 W ab + t(t b t a ) JZ H 2 W ab + ^(L^l^G^ 

k=i 



(7.43) 
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The color structures Ck are denned in Eq. (17.61) and there are six new structures (£fc)aL : 
L\ = (t a t b t c )ji, L 2 = (t c t a t b )ji, L 3 = (t a t c t b )ji, 

U = {tHh 11 )^ , L 5 = (tH^)^ , U = {t b t c t a ) 3i . (7.44) 

The invariant kernels for the first operator, X l j B = ip^ip^ ', are given by 

[Hlflizt, Z 2 ) = ~ Z nlj da [ d PP *l2i Z 2l) 

+ / da I dp I dmp(z" 2 , z% % , z^) , . 



o Jo J p 

1 pot pi 



[H 2 ip](zi, z 2 ) = -z 12 I da j dp I ^77 <p(z" 2 , z p 2X , z^) 



[H 3 (p] (z u z 2 ) =^2<J da J df3^(p{zi, z 21 ,z( 2 ) 



da ["dp f d^ip(z^z^zl) } . 17. 15) 

Jo Jo 



»1 pi 

[HM(zi,z 2 ) = z\ 2 I da dPa 2 (3cp(z" 2l z^) , 

JO Ja 

[HM(z u z 2 ) = z\ 2 J da £ dPaap (2 - ^ <p(z 21 , 4) • (7-46) 



and 



[Hnp]{zx,z 2 ) = z\ 2 I da j dpi c?7 ^{z^z^z^) 

JO Ja J P 



a 

[H 2 <p}(zi, z 2 ) = z\ 2 da j dp apj \2 - <p(z? 2 , z^, z^j , 



[Hzip\{z u z 2 ) = z\ 2 I da! dp! dj (p(z? 2 , z^, zl x ) 



J a J a 

-a p/3 



[H 4V ](z u z 2 )=z 3 12 J^ da£dpj o ^^7(1 + ^(1-^)) <p(*&,4i, 

pi pa. pi 

[H 5 (p](zi,Z2) = z\ 2 I da j dpi dj ap^f ip(z X2 , z^, z^) , 

Jo Jo Ja 

[H 6 <f}(z 1 ,z 2 )=zf 2 [ da! dp! djaPj^ 2 ,z^,zl). (7.47) 

Jo Jo J/3 
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For the second operator, X l j[ B = ^D^ + ip+ A ip 3 + B , we obtain 
[Hip] {zi,z 2 ) = z 12 



1 f 1 f 1 P 2 

dp P (p(z u z 2 , 4) + / da I d(3 — (p(z 21 , z 2 , z{ 2 ) 



JO J a 

+ 2 / da / d/3 / d7 7#i tt 2 ,4>4) 

'0 JO J/3 



[W 2V 7](^, a*) = zj j 1 da J X d(3^ (^2 - ^ 22 ^ 



5 ^21, 



/*1 /*tl /»1 

+ 2 / da / d/3 / ^7^(4-4,4) 
Jo Jo 

-1 /-a 



[H 3 <p}(z u z 2 )=2z 12 [ da f dp f dyytp(z? 2 , 4, 4) , (7.48) 

JO JO JO 

[Hi(p](zi,Z2) = -Z12 J da J dfiaafi ^2 + (p(z^ 2 , z^) , 

[H 2 <f}( Zl ,z 2 ) = -z 12 J^da j_d(3aaP ( 4 ~ ) ^i,4) ^ 49 ) 



and 



[Ki<p](zi, to) = -z$2 da J dp J dj aPy \2 + ip(z^ 2 , z^, z^) ■, 
[H 2 <f](zi, z 2 ) = -z\ 2 1^ da j dp dy apy ^4 - <p(z" 2 , 4, 4) 



[H 3 <p](zi,z 2 ) = -zf 2 da J dp J dj aPj (2 + ip(z^ 2 , z^, z^) , 
\n^\{z^z 2 ) = -z\ 2 f da f dp [ d 1 aP 1 (2 + ^(2-^X\ y?(^ 2 ,4,4) 



JO JO 

1 pa r l 



da J dp J dy apy ( 1 + + -=-) <p(z? 2 , 4, 4) , 
[H 6 <p]{zi,z 2 ) = -zl 2 ^daJ\/3£dya C P'y (l + °t + t^j ^2,4,4)- (7-50) 

8 Conclusions 

Extending the work by Bukhvostov, Frolov, Lipatov and Kuraev (BFLK) |32j, in this 
paper we derive a complete set of two-particle 2 — > 2 RG kernels for collinear twist E = 3 
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operators and also the 2 — > 3 mixing kernels into three-particle quasipartonic operators, 
cf. Fig. 1. The kernels are written for the renormalization of light-ray operators built of 
chiral fields in a particular basis such that the conformal symmetry is manifest. They 
serve as building blocks in the RG equations for arbitrary gauge-invariant twist-four 
operators in Quantum Chromodynamics. The results can easily be recast in momentum 
fraction space, in the form of evolution equations for higher-twist generalized parton 
distributions (cf. [3 HE]). Specific applications will be considered elsewhere. 

To this end we suggest a new technique which is based on using mostly algebraic 
methods and bypasses calculation of Feynman diagrams. The main idea is that RG 
equations for operators containing field components of different collinear twist are re- 
lated by Poincare symmetry. In this way the SX(2, R) invariance of BFLK kernels for 
quasipartonic operators becomes extended to the full conformal group SO (4, 2). As the 
result, the 2 — > 2 RG kernels involving "minus" quark and/or gluon fields can be ob- 
tained from the results of Ref. [32J by a simple replacement of the two-particle quadratic 
SX(2, R) Casimir operator by the SO(4,2) one, cf. Ref. [37]. The way to explicit ex- 
pressions for 2 — > 3 kernels is somewhat longer: They can be obtained by applying 
translations and Lorentz rotations in the transverse plane to the renormalized leading 
twist operators. We found this technique to be very efficient. It also offers a new insight 
in the algebraic structure of operator renormalization in gauge theories and may be use- 
ful in a more general context. It would be interesting to investigate whether the same 
methods can be used beyond the one-loop approximation. 
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Appendices 



A The Conformal Group 

Generators of the conformal group act on the fundamental fields in the spinor represen- 
tation in the following way [49J: 

i[P a6c ,$(x] 

' M,, ;• 'H- r 

i[Kaa, ${x] 



d aa <$>(x) = iPaaHx) , 

~ (xaiff* + 2t + r ^ + Cdi) *(*) = %D $ (x) , 

| (x^di> + x^aj - 2£ a A - 2^^) $(z) = ^$(x) , (A.l) 



r 9 



2£c 



d£f> 



9 



$(x) = iM d 3$(x) 



9 



Here $ = $ e ) with $ s = {V'frXfJ&h = {V'C.Xfc/tf} where £ is an auxiliary 
spinor. The generators P^, P, M a/3 (M d «) and P"^ correspond to translations, dilata- 
tion, Lorentz rotations and special conformal transformations, respectively. 

The chiral fields transform according to the (s, 0) representation of the Lorentz 
group 



2^0) 



1)$ 



(SO) 







2 - - c ' , 2 ^ 

where M 2 = M Q ,.\/'"'. M 2 = .\ /. ,.!/" ' and s = 1/2 or s = 1 for quarks and gluons. In 
turn, the antichiral fields transform according to the (0, s) representation 



-ilP$ 



2^05) 







1 ™ 



-ilP$ 



2* (OS) 



s(s+ 1)$ 



(OS) 



2 5 ' 2 

Finally i = 1 is the geometric twist [JT] of the field $: It is defined as t = £ can — s — s 
where £ can is the canonical dimension. 

The commutation relations for the generators read 

[pad D ] = _ ipa d ^ ^od D ] = ^ad ? 



[-Pad) -^87] 
[M^ M jS ] 



[^a/^-Pya ^a"fPf3a\ > [-fad) 



/3 7 J 



[ e o/3-^a7 + e d7^a/3] ' 



(e 7/3 M a5 + e 5/3 M a7 + e ia Mp S + e^M^) 



2 ( e /?7^d«j + tyMv, + eajM $ s + e^M .) 
4* (siSiD + siMj' + W 



(A.2) 
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The expression for the quadratic 50(4, 2) Casimir operator C 2 reads: 

C 2 =i {K a «P a « + P a6l K a ^ ~\d 2 + 1 - (m^M^ + M^M* ) ■ (A. 3) 

One finds 

C 2 $( s0) = \(s 2 - 1) $^ 0) , C 2 $f 5 = ^(s 2 - 1) <lf } . (A.4) 
The generators of the collinear SX(2,IR) subgroup are defined as 
S + = % -{jiKp), 5_ = ~(APA), 

5o = ^ ((M)(A»£ - (Xfi^M^ - (fi\)M & ^\^ , (A.5) 
In order to shorten the following expressions, without loss of generality one can assume 

m\ 



X a = (1,0), X a = (0,1), 

//* = (0,1), /i Q = (-l,0), (A.6) 
so that (/zA) = (A/i) = 1. Using this convention the same generators can be written as 

S + = \K 22 , S_ = -~P xi , S = % - [D - M 12 - M i2 ) . (A.7) 

The SX(2,K) generators commute with the operators of the collinear twist E and chi- 
rality H defined as 

E = i(D + M 12 + M i2 ) , H = i (M i2 - M 12 ) . (A.8) 

For our purposes it is useful to separate the contribution of the light-cone subgroup to 
the expression for the two-particle Casimir operator C 2 2 : 

C 2 2 - C 2 - C 2 = (S x + S 2 ) 2 -S 2 -S 2 2 + X -E X E 2 + H X H 2 

+ \ (Kl 2 P 2>12 + K 2i P 2j21 + K l 2 2 P hl2 + K 2i P li2i + K 22 P 2i22 + K 22 P h22 ) 

+ (M 22 M 2 n + M\ l + M 22 M 2 1{ + M^M 22 ) , (A.9) 

where e.g. P 2 12 stands for the (12) component of the translation operator acting on 
the second field, etc. The terms in the last two lines all vanish if applied to the tensor 
product of two "plus" light-ray fields, so that one finds 

C\ 2 ® = (S 2 12 - k 12 ) ® $f . (A.10) 
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Here Sf 2 = (Si + S2) 2 is the two-particle Casimir operator of the collinear SL(2, 
subgroup and K± 2 is a constant: 



«12 



2 - \ (hi + h 2 f 



(A.ll) 



where /ii, /12 are the helicities of the fields: H k $ 



(*) 



B Invariant kernels 

In this appendix we formulate a simple rule how to construct SL(2, M.) invariant kernels. 

First, note that if a function ^(^1,22) transforms according to the representation 
T jl <g> T j2 of the SL(2, R) group, i.e. 

, 1 /azi + & az 2 + &\ m 1 x 

then the function f(zi,z 2 ) = (z\ — z 2 ) 2n ip(zi, z 2 ), n being integer or half-integer, trans- 
forms according to the representation T J1_n ® T- 72- ™. An invariant operator 7i which 
maps T jl <g) T- 72 — > T n ® T 12 exists if and only if ji + j 2 — *i + ^2 + 2m, where m is, 
again, integer or half- integer. It is clear that invariant kernels for m/fl, H,( m \ can be 
represented in the form Ti.^ = (z\ — z 2 )' 2m l-L^ where TC is an invariant kernel for the 
case ji + J2 = H + It is, therefore, sufficient to consider this case only. 

To write down the most general expression for an invariant operator 7i : T- 71 ® T- 72 t— > 
Jin ^ j"i2 ac ting on the space of functions of two variables to 2 ) l— ► vK^i) ^2)1 consider 
the following diagram: In this diagram the points ti>i,ii>2 are connected with zi,z 2 by 




"22,^22 



four lines. Each line, e.g. connecting the points Z{ and Wk, carries a pair of indices, 
(a ik ,p ik ), i, k = 1, 2. The indices on different lines are subject to the constraints: 

otik + a 2k = l, p lk + p 2k = 2j k , p kl +p k2 = 2i k , k = l,2. (B.2) 

An invariant kernel has the form [17] 

[H(f](zi, z 2 ) = da 21 da l2 T\ a^ 1 u f a2l0!l2 ] <p( Wl (a), w 2 (a)) , (B.3) 
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where 

2 

m=l 

and u(x) is an arbitrary function. 

The choice of p-parameters which satisfy Eq. ( IB. 21) on is not unique, but the ratio of 
the integration weights corresponding to two different solutions is a power of the invariant 
(anharmonic) ratio (a 2 i(Xi 2 ) / (ceua 22 ) , i- e - a different choice results in the redefinition of 
the function to. 

Changing the notation to 

an = 1 - a , a 21 = a, a 12 = (3, a 22 = 1 - /? , Pn = «, P21 = a, P22 = b, p u = b 
Eq. ( IB. 31) takes a more familiar form 

[Wd(*i,2a) = J^a J^dPa*- 1 ^- 1 ^^^ {^j tp{*&,2&) ■ (B.4) 

Let us compare this general expression with the kernels given in (14.11) . Notice that the 
integrations in (14.11) do not go over the unit square, < a, (3 < 1 as in (1B.4I) . but over the 
simplex either 0<a + /5<lorl<a; + /3. These integration regions are not arbitrary, 
since the only way to obtain them is to choose the function u(x) proportional to the 
theta-function of the anharmonic ratio. In particular, the choice u(x) = 8(1 — x)ui(x) 
results in the integration region < a + (5 < 1. The kernels in (14.11) that contain one- 
dimensional integrals arise when w(x) is proportional to a 5— function. For example, 7i d 
corresponds to the choice u(x) = 5(1 — x). The 'exchange' kernels correspond to 

the solutions with 

a = 2ji - k - 1 , a — k + 1, b = 2j 2 - 1 , b=l 

and u(x) = 5(x). 

More complicated, e.g. 2^3 kernels H : T jl ® T j2 ® T js ^ T h ® T h , can be 
constructed similarly. It is again sufficient to consider the case ji + j 2 + J3 = i\ + i 2 - To 
this end, draw a diagram with wi,w 2 ,W3 points connected to each of the Zi,z 2 points: 
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Similar to the 2 — > 2 case we furnish each line with a pair of indices which have to 
satisfy the following equations: 

aifc + Qf2fc = l, fc = l,2,3 

Pik+P2k = 2jk, A; = 1,2,3 

Pml + Pm2 + Prn3 = , m = 1, 2 . 

The most general invariant kernel 7i has the form 

3 pi 2 3 / ^ \ 

[H<p]{z 1 ,z 2 ) = TT / da lk TT TT of™" 1 u( 21 12 , 21 13 ) <p(w 1 (a),w 2 (a),w 3 (a)) , 

H JO f=lm=i \ a ll a 22 ail«23/ 

where Wfe(a) = Ylm=i z m®-mk and a;(x, y) is an arbitrary function of two variables. 
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